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The effect of topography on the evolution of a disturbance in a baro- 
clinically unstable mean flow is studied using analytical and numerical 
Simulations. In particular, dynamical mechanisms involved in cyclogene- 
sis in the lee of long, meridional barriers similar to the Rocky Moun- 
tains are explored. The rapid growth observed in lee cyclogenesis is 
highly dependent on the superposition of a growing baroclinic wave with a 
steady, orographically forced wave of the same scale. As the baroclinic 
wave moves over the mountain, development is masked by the orographical- 
ly-forced, high-pressure ridge. As it moves down the lee side of the 
mountain, the baroclinic wave appears to grow rapidly due to superposi- 
tion with the forced, lee-side trough. Indications of enhanced, lee-side 
baroclinic instability are present in the numerical simulations, but the 
effect on wave development is minor compared to the effect of superposi- 
ae The rapid continuous-mode growth which has recently been demon- 


strated by Farrell (1982) is not observed in any of the simulations. 
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Steady, forced solution Y =Y m(z) cos ux given by (2.13). 27 
Solution is given at: z = 0.0 (solid curve), 

z = 0.5 (dotted-dashed curve), and z = 0.1 (dashed 

curve); z = - 1n p/ps is the non-dimensional 

vertical coordinate. The shaded region at the 

bottom represents the terrain, h = hs cos UX. 


Amplitude of the steady forced solution v = 3px 27 
as a function of non-dimensional height, z, for 

global wave numbers: uy, = 8 (solid), u, = 12 

(dotted-dashed), and yw, = 16 (dashed). 


Steady, forced solution for the linear model presented 29 
as Pg, the deviation from mean-state surface 

pressure, for a 45% sector if the globe. Cyclic 

continuity is assumed at the east and west boundaries 

of the sector. Profile of the mountain given by 

(2.14) is included at bottom. 


Transient, free solution for the linear model given 30 
by (2.1, 2.5) presented as pg for a 450 sector of 

the globe. Curves are given at six-hourly time 

intervals at 450N. Cyclic continuity is assumed at 

the east and west boundaries of the sector. 


As in Figure 2.4, except for combined (free + forced) 31 
linear model solution. Profile of the mountain 
given by (2.14) is included at the bottom. 


Minimum pe versus time at 45°N for free (solid) 33 
and combined (dashed) linear model solutions. 


Phase versus time for free (solid) and combined 34 
(dashed) linear model solutions. Phase is inferred 

an the location of the minimum value of pg at 

459N. 


Steady, forced solution of the semi-geostrophic model, 45 
VM» in transform space. At the points -X1 and 

+X1, 9Vm/9X = 0; these points separate the region ` 

of phase acceleration from the regions of deceleration. 
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Semi-geostrophic model solutions presented as p< 
along 450N in transform space (solid) and in 
physical space (dashed). The 3.0km mountain is 
included at tne bottom. 


Horizontal distribution of variables in the UCLA 
model grid. 


Vertical distribution of variables in the six- 
layer version of the sigma-coordinate system. 


Meridional cross-sections of the initial mean wind, 
u, for: a) 160 jet, and b) 80 jet. Contour 
interval is 10 m s”?. 


Plan-view of the topography on a 45% sector of the 
globe. Vertical line at the center of the figure 
represents the mountain ridge line. 


Amplitude of wave number 8 of v at 45%N as a function 
of yg for linear experiment 1. Initial amplitude 
(dashed) versus o is included. 


Amplitude of wave number 8 of v ato = 0.9167 at 
45°N plotted as a function of time for linear 
experiment 1. 


Amplitude of wave number 8 of v at 45%N as a function 
of o for linear experiment 2. Initial amplitude 
(dashed) versus o is included. 


Amplitude of wave number 8 of v at 45°N as a function 
of o for linear experiment 3. Initial amplitude 
(dashed) versus o is included. 


Vertical cross-sections of v along 45°N at 12 h from 
the six-layer (solid) and 20-layer (dashed) models 
for linear experiments 1, 2, and 3 (a, b, and Es 
respectively). Contour interval is 2 m s”?. 


Sea-level pressure contours for nonlinear Experiment 1 
Aid eae lS nec) 24h. d) 30h, e) 36 h, 

f) 42 h, and g) 48 h. Contour interval is 2.5 mb. 
Domain is a 45% sector of the globe with cyclic 
continuity at the east and west boundaries. 
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Amplitude of wave number 8 at g = 0.9167 at A5ON 
versus time for nonlinear Experiments 1 (solid) 
and 2 (dashed). 


Amplitude of wave number 8 at o = 0.9167 at 45°N 
versus time for nonlinear Experiments 3 (solid) 
and 4 (dashed). 


Sea- level pressure contours for control run A at: 
a) O hy bi 24 nee) 4S ande SDE 
Contour interval is 2.5 mb. 


Amplitude of wave number 8 of v at o = 0.9167 
versus time at 45°N for control run A. 


Amplitude of wave number 8 of v ato = 0.9167 
versus time at 450N for control run B. 


Sea-level pressure contours for the Interaction 
and Non-interaction solutions of Experiment I, 
Case 1. Interaction solutions are given at: 

a) to» bD) totl2 h, c) tpy+24 h. The non- 
interaction solution dl is given at tot24 h. 
Contour interval is 2.5 mb. Bold contours at 
center of figure represent mountain. 


sea-level pressure contours at t¿+24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment I, Case 4. Contour interval is 

2.5 mb. Bold contours at center of figure 
represent mountain. 


Amplitude of wave number 8 of v at o = 0.9167 at 
45°N for Experiment I, Cases 1, 2, 3, and 4. 
Interaction (dashed) and Non-interaction (solid) 
solutions are shown. 


Minimum pé (where pc = Ps-P<) at 450N versus time 
for Experiment I, Cases 7. 3° 3, and 4. Interaction 
(dashed) and Non-interaction (solid) solutions are 
shown. 


p< (where ps = Ps-Ps) versus X for Experiment I, 
Case 1 at 45°N at 6 h intervals from ty to t +24 h. 
Dashed curves represent Interaction solutions; 
solid curves represent Non-interaction solutions. 
East-west profile of mountain is shown at bottom. 
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Ps (where ps = ps-ps) versus X for Experiment I, 
Cases 2, 3, and À at 450N at 6 h intervals from 

to to tot24 h. Dashed curves represent Interaction 
solutions; solid curves represent Non-interaction 
solutions. East-west profile of mountain is shown 
at bottom. 


Sea-level pressure contours at t)+24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment II, Case 1. Contour interval is 
2.5 mb. Bold contours at center of figure 
represent mountain. 


Amplitude of wave number 8 of v at go = 0.9167 at 
450N for Experiment II, Cases 1, 2, 3, and 4. 
Interaction (dashed) and Non-interaction (solid) 
solutions are shown. 


Sea-level pressure contours at tpot+24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment III, Case 1. Contour interval is 
2.5 mb. Bold contours at center of figure 
represent mountain. 


Sea-level pressure contours at t +24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment III, Case 4. Contour interval is 
2.9 mb. Bold contours at center of figure 
represent mountain. 


Amplitude of wave number 8 of v at o = 0.9167 at 
450N for Experiment III, Cases 1, 2, 3, and 4. 
Interaction (dashed) and Non-interaction (solid) 
solutions are shown. 


Sea-level pressure contours at t)+24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment IV, Case 1. Contour interval is 
2.5 mb. Bold contours at center of figure 
represent mountain. 


Sea-level pressure contours at tp)+24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment IV, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure 
represent mountain. 
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Sea-level pressure contours at to+24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment Va, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure 
represent mountain. 


Position of surface-pressure trough at to, 
to+12 h, and to+24 h in Experiment Va. 


Sea-level pressure contours at t +24 h for 
Interaction (a) and Non-interaction (b) solutions 
of Experiment Vb, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure 
represent mountain. 


Position of secondary surface-pressure low at 
three-hourly intervals for Experiment Vb. 
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I. INTRODUCTION 


In recent years, there has been renewed interest in the effects of 
topography on atmospheric flows. Mountains of various shapes, sizes and 
orientations force motions on all scales - from small, turbulent eddies 
to planetary waves. The objective of this research is to determine the 
crucial mechanisms in lee cyclogenesis, wherein cyclones form on the 
downwind side of major mountain ranges. 

It has been generally accepted for a number of years that a dynamical 
connection exists between the presence of mountains and the occurrence of 
lee cyclogenesis. This is based on evidence presented in both observa- 
tional and numerical studies. Petterssen's (1956) climatological statis- 
tics of cyclogenesis frequency established that the lee sides of large 
mountain complexes such as the Rockies, the Alps and the Himalayas are 
preferred regions for cyclogenesis. More recently, the observational 
Studies by Reitan (1974); Chung, Hage and Reinelt (1976); Zishka and 
Smith (1980); as well as others, have produced similar findings. In 
addition, various numerical studies, such as those by Egger (1974); 
Manabe and Terpstra (1974); Bleck (1977); and Tibaldi, Buzzi and Malguzzi 
(1980), have concluded that inclusion of topography in numerical models 
is necessary to predict cyclogenesis in the "preferred" regions. 
Consequently, it is generally accepted that mountains play a role in the 
cyclogenesis process. However, the dynamical nature of that role is 


still not well understood. 


PS 


The basic components of lee cyclogenesis appear to include: large 
mountains, the presence of the jet stream normal to the mountain barrier, 
and a pre-existing disturbance. The required mountain size is somewhat 
arbitrarily defined, although the results of Petterssen (1956) and others 
indicate that cyclogenesis is not observed with enhanced frequency in the 
lee of mountains with smaller heights and widths, such as the Urals. The 
presence of the jet stream during lee cyclogenesis has been thoroughly 
documented in a number of studies, including Newton (1956), Klein (1957), 
Hovanec and Horn (1975), and Whittaker and Horn (1980) for the Rockies, 
and Trevisan (1976) and Buzzi and Tibaldi (1980) for the Alps. In the 
typical case (see, for example, Hess and Wagner, 1945; Palmen and Newton, 
1969; and Buzzi and Tibaldi, 1980), lee cyclogenesis is initiated as a 
disturbance crosses a large mountain range. During the initial stages of 
development, rapid intensification is observed with e-folding times of 12 
hours or less (Buzzi and Tibaldi, 1980). This e-folding time is much 
less than the rough estimate of 25-30 hours given for cyclogenesis 
occurring away from mountains (Palmen and Newton, 1969; Holton, 1971). 
After moving away from the mountain, the growth and phase speed of the 
lee cyclone approach those values normally associated with cyclone de- 
velopment. Consequently, it has been suggested that topography alters 
the baroclinic instability process in a way that results in more rapid 
intensification on the lee side of large mountain ranges. 

The basic objective of this research is to gain a better understand- 
ing of how topography affects the cyclogenesis process. This investiga- 
tion will provide insight into the dynamical basis for the rapid 


intensification which is observed on the lee side of large mountain 
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ranges, and, in so doing, provide a better understanding of why these 
areas are preferred regions for cyclogenesis. Since the nature of the 
orographic effect is dependent on the size, scale and orientation of the 
mountains (Palmen and Newton, 1969; Chung et al., 1976; and others) the 
types of mountain ranges to be considered must be limited. This research 
focuses on cyclogenesis in the lee of a fone meridional barrier which is 
intended to simulate the Rocky Mountains. Cyclogenesis is treated as a 
manifestation of ie instability, and orographic effects are 
sought which can explain the increased growth which is observed on the 
lee side of large mountains. 

Recent dynamical investigations of flow over topography provide a 
number of possible mechanisms which could play a role in the cyclogenesis 
process. Charney ‘and DeVore (1978), and Charney and Strauss (1979) con- 
Sidered topographic forcing of atmospheric flow and identified the form- 
drag instability which results from the interaction of mean flows with 
large-scale mountains. This mechanism, which depends on the Rossby wave 
being nearly stationary, appears to play a role in the evolution of 
planetary waves and blocking patterns; these waves occur on scales larger 
than those observed in lee cyclogenesis. Pedlosky (1981) considered 
resonant topographic forcing as a growth mechanism and found that slowly 
moving long waves could experience large increases in amplitude due to 
resonance. This mechanism also appears to be applicable to the longer 
and more slowly moving planetary waves rather than the cyclone waves, and 
it is not investigated here. 

Merkine (1975, 1977) and Merkine and Israeli (1978) have considered 


the interaction of baroclinic flows with topography using the 
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quasi-geostrophic and semi-geostrophic equations. Merkine's (1975) 
solution for a vertically sheared, steady, quasi-geostrophic flow past an 
infinitely long ridge showed that the flow acquires a vertically sheared 
northerly component of velocity along the lee side of the mountain. This 
vertically sheared northerly component is in addition to the shear of the 
basic state and indicates enhanced baroclinicity. According to Merkine, 
the increased baroclinicity could explain the higher growth rates 
observed in os However, it has not been demonstrated how 
such enhancement could increase the instability of lee cyclones. 

Farrell (1982) used the linearized quasi-geostrophic equations to 
demonstrate that continuous-mode solutions can, with the proper initial 
conditions, exhibit large growth rates for periods of up to 12 hours 
after excitation; this initial growth is followed by damping. The addi- 
tion of these modes to the exponentially growing discrete modes could 
explain the enhanced growth which is observed in lee cyclogenesis. The 
selection of initial conditions is crucial for this effect, since one 
could always pick a special initial state with will give rapid initial 
growth. 

Smith (1984) has proposed a baroclinic lee wave mechanism for lee 
cyclogenesis in the Alps. This mechanism depends on a mean flow in which 
the vertical wind shear is opposite to the surface component across the 
mountain range. However, it does not appear to apply to the Rocky Moun- 
tain cases where the vertical shear is usually in the same general 
direction as the surface wind. 

Still another possibility is that the mountains have only an apparent 


effect on the cyclogenesis process. A steady flow over a synoptic-scale 
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mountain forces a stationary high-pressure ridge over the mountain and a 
pressure trough on the leeward side (see, for example, Smith, 1979). A 
growing cyclonic disturbance of approximately the same scale as the moun- 
tain would be cancelled by superposition with the high pressure ridge as 
it moves over the topography. On the lee side, large growth rates would 
be observed as the cyclonic disturbance continues to grow and become 
superposed with the lee-side trough. With a high enough mountain, such 
lee-side growth could appear as rapid deepening. In addition, the lee- 
Side superposition would induce a closed circulation at an earlier stage 
of the amplification process than would occur over the mountain or flat 
terrain. As a result, the initial appearance of the closed surface- 
pressure contour (which is defined as cyclogenesis in most of the studies 
quoted above) would be accelerated by the presence of the mountain. This 
could explain the higher frequencies of lee cyclogenesis, since growing 
waves placed randomly are more likely to have the first closed contour 
over the lee than over the mountain top. 

Of the dynamical mechanisms described above for lee cyclogenesis, the 
three most likely are examined in this study: 1) enhanced, lee-side baro- 
clinic instability; 2) continuous-mode growth; and 3) superposition. 
Chapter II treats analytically the interaction between an infinitely long 
mountain range and baroclinic waves in the basic state which was used by 
Eady (1949). The quasi-geostrophic equations for this problem are linear 
So that the mountain-forced solution and the unstable Eady solution are 
mathematically independent. An example is given which shows that the 
Superposition of the two simple solutions can give a reasonable represen- 


tation of lee cyclogenesis. Semi-geostrophic solutions are also 
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derived which show some interaction between the mountain and the moving 
disturbances. This interaction changes the phase speed of the moving 
wave, but it does not change the intensity of the disturbance. 

In Chapter III, the finite-difference and spectral numerical models 
which will be used in later chapters are described. The topography and 
initial conditions for the following numerical experiments are found in 
Chapter IV. In Chapter V, the continuous-mode growth mechanism which was 
proposed by Farrell (1982) is investigated. The initial conditions are 
crucial to this mechanism. In this study, simple orographic forcing is 
used to perturb a basic state that initially contains no perturbations. 
Such a situation might be expected if the mean flow was rapidly increased 
from zero to the limiting value. The linear study in this chapter uses 
the spectral model and the nonlinear study uses the finite-difference 
model. None of these experiments show any significant continuous- 
spectrum growth. 

Chapter VI contains the numerical experiments with more realistic 
conditions which include a pre-existing disturbance and a baroclinic jet 
with vertical and horizontal shear. A procedure is developed which 
allows a comparison between a transient mountain solution and a solution 
in which a no-mountain transient solution is superposed on a forced moun- 
tain solution. The experiments are carried out for a range of initial 
disturbance amplitudes so that the small-amplitude experiments show the 
linear behavior of a wave in a mean flow forced by the mountain. The 
larger amplitude experiments include nonlinear interaction between the 
disturbance and the mountain-forced flow. The experiments show no signi- 


ficant difference in intensity between the mountain interaction 
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cases and the superposition cases. Phase changes are found in the inter- 
action cases which are consistent with those predicted by the semi- 
geostrophic theory in Chapter II. The various studies included in this 
research lead to the conclusion that superposition is the main mechanism 
for lee cyclogenesis with mountain systems similar to the Rocky 


Mountains. 
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IT. ANALYTICAL MODELS 


Baroclinic instability has been accepted as a major source of propa- 
gating synoptic-scale disturbances in the mid-latitudes. The success of 
this theory in predicting the general scale, growth, and structure of 
cyclone waves is undisputed. In this chapter, some simple analytical 
models are introduced to provide a theoretical basis for investigating 
the dynamical mechanisms involved in lee cyclogenesis. 

In the first section, the model developed by Eady (1949) is briefly 
reviewed. Although this model does not include topography, it provides a 
good description of wave development in the early stages of cyclogenesis 
and is useful for examining the linear dynamics of baroclinic instabil- 
ity. In the second section, Eady's model is modified to include small- 
amplitude topography through the lower boundary condition; the forced, 
Steady solution is derived for investigating orographic effects on wave 
development. In the third section, a non-geostrophic model is used to 
describe wave development in the presence of a finite-amplitude mountain. 
The potential vorticity equation for this model is transformed following 
Hoskins and Bretherton (1972) and solved analytically. Analysis of the 
solution reveals that finite-amplitude topography distorts the structure 
of a developing baroclinic wave while in the vicinity of the mountain, 
but has no net effect on growth away from the mountain. As will be 
demonstrated, these simple models appear to explain very well some of the 


dynamical features observed in lee cyclogenesis. 
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A. REVIEW OF THE EADY MODEL 

Eady (1949) considered a simple, yet realistic baroclinic instability 
problem. His model uses the linearized form of the quasi-geostrophic 
equations with the Boussinesq approximation to describe an inviscid atmo- 
Sphere confined between two rigid horizontal plates. The basic problem 
is stated in terms of a single partial differential equation which 
expresses conservation of potential vorticity. The vertical velocity is 
assumed to vanish at the boundaries, z = 0, 1 (where z = - 1n p/po is the 
non-dimensional vertical coordinate). The following assumptions are 
made: The basic wind flow varies linearly with height - that is, U = Sz, 
where S is a constant. The supporting horizontal temperature gradient is 
constant. The basic state density, ©, and the static stability, fr, are 
taken as constants, and the effect of the earth's sphericity is ignored 
by assuming that the Coriolis parameter is constant. With these assump- 
tions and the normal mode solution form 
e -ct 2 Malo) 


Y > Re 


the quasi-geostrophic potential vorticity equation reduces to the ordin- 
ary differential equation for Y¥r(z) given by Haltiner and Williams 


(1980): 


d? 
E oar Fee cc! = 0 
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The subscript, F, is used here to denote the free, or unforced, solution. 
The boundary conditions are obtained from the first law of thermodynamics 


d o, əp aw dU 


weet al 2 so (203a) 
O 


(r y u. dz  ðx əz 


If (2.1) is substituted into (2.3a) and the non-dimensional vertical 





velocity 
2 = Re [weze Oct , then (2.3a) becomes 
(U N c) otg = Uy + TW = 0 i (2.3b) 
82 IS 


At the upper and lower boundaries where Z = 0, the boundary conditions 


become 
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The free solution takes the form 
Ye = A sinh (£) +B cosh (£) (2.5) 
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Application of the boundary conditions yields the characteristic equation 
for the phase speed, c. The condition for baroclinic instability is that 


c must be complex, which occurs if 


and 570 : (2.6) 





where L = 2 7/1. When these conditions are satisfied, Y in (2.1) becomes 
an exponential function of time, with growth rate proportional to verti- 


cal shear. Maximum growth rate is predicted for 





this indicates that wavelengths of the order of 4,000 km will dominate 
the spectrum of atmospheric waves, and is in approximate agreement with 


the observed scale of mid-latitude cyclones. 


B. EADY'S MODEL WITH SIMPLE TOPOGRAPHY 

The effect of topography can easily be included in Eady's model if 
the topography has a small amplitude compared to the scale height of the 
atmosphere. In this case, the basic flow remains zonal and the effect of 
the topography is felt through the lower boundary only. This simplifica- 
tion is frequently made for mathematical convenience and provides a rough 
approximation of the effect of the topography on the flow. As will be 
Shown, even with this weakness, the model is able to describe at least 
qualitatively the increased growth observed on the lee side of the 


mountain. 
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1. Model Formulation 
Consider the case of steady flow over a small-amplitude mountain 


with a profile given by 


h(x) = (re 1,00) peo 


where he is a constant. Notice that h is independent of y, which indi- 
cates a mountain of infinite north-south extent. Let the basic flow be 
given by 


U = Us + Sz me) 


where Us, the surface wind, is constant. It will be shown below that the 
inclusion of the surface current is necessary to obtain an orographic 
effect through the lower boundary. The surface current has no effect on 
the preferred scale and growth rate of the basic Eady solution reviewed 


above. If the solution form 


2 ju(x - ct) $ (2.9) 
y = Re (Ye 
is substituted into the steady form of the quasi-geostrophic potential 


vorticity equation, the equation becomes 


fo d “ty 
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where the subscript, M, indicates the steady, mountain-forced solution. 
The boundary conditions are obtained from (2.3b). At the upper boundary, 


where z = Q, 


dy , (2.11a) 
(Us + S) 2 - Sy =0 
while at the lower boundary, 
dY h 
A EA > (2.11b) 
S dz C 


where H is the scale height of the atmosphere. The term in (2.11b) which 
includes hs follows from the topographically-forced vertical velocity at 
the lower boundary 

ðh À (2.12) 


. S 
Z(x,0) = — — 
| ) H 9x 


The solution of (2.10) with boundary conditions (2.lla, b) is 


E 
w Z Z 1 
SJ cos (2.13) 
L <2 2 
F H(U a 2 y) € € 
where 
f? 
O -33 
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This result shows that the amplitude of the forced wave is proportional 
to the magnitude of the surface wind and the height of the mountain. As 
indicated previously, a non-zero surface wind is required for an oro- 
graphic effect on the flow. The solution yp = Yue UX is shown for various 
levels in Figure 2.1. A high-pressure ridge is centered over the moun- 
tain with troughs on both the upstream and downstream sides. The phase 
is constant with height, but the amplitude damps exponentially. The 
amplitude of the meridional wind component, v = 9W/dx is shown versus 
height for various wave numbers in Figure 2.2. The amplitude of v damps 
in the vertical in each case, and the magnitude of the damping increases 
with wave number. This latter result indicates that, as the horizontal 
scale of the mountain is decreased, the vertical extent of the orographic 
effect is reduced. These solutions are generally similar to those given 
by Smith (1979) for an atmosphere with no vertical wind shear. 
2. Model Application 

The steady, forced solution given by (2.9) may be combined with 
the free Eady solution to form a linear model of a wave passing over 
topography. As a result of the linearity assumption, the free and forced 
solutions do not interact. Therefore, this model does not predict any 
new or enhanced instability mechanism associated with the mountain. How- 
ever, due to superposition, a developing, small-amplitude, baroclinic 
wave of the proper scale can be cancelled by the forced ridge as it moves 
over the crest of the mountain and appear to grow rapidly as it moves 


down the lee slope and into phase with the stationary trough. 
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Figure 2.1. 


igure 2.2. 





Steady, forced solution Y = Y m(z) cos ux given by (2.13). 
Solution is given at: z = 0.0 (solid curve), z = 0.5 
(dotted-dashed curve), and z = 0.1 (dashed curve); 

Z = - In p/p. is the non-dimensional vertical coordinate. 
The shaded region at the bottom represents the terrain, 

h = hg COS UX. 


Amplitude of the steady forced solution v = dW/dx as a 
function of non-dimensional height, z, for global wave 


numbers: Uy = 8 (solid), uy = 12 (dotted-dashed), 
and u, = 16 (dashed). 
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Let the constants in equations (2.1 - 2.13) be given the follow- 
ing mid-latitude values: fp = 10-4 s-l, r = 10%més-¢, Us = 5 m s“l, and 
S = 40 m s-1. The constant B in (2.5) represents the initial amplitude 
of the free wave and is arbitrarily assigned the value of 300 mês Let 


the mountain profile be given by the truncated Fourier series: 
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z ; 2.14 
h(x) = ZA = ) A cos nx) a) 
n=1 
where 
Ag = 1/4 (2.14b) 
ASIA 2 (2.14c) 


A 
Nn  nw(4-n?/4 


and z, = 1.5 km. Ao is evaluated by taking the limit of Ap as n 
approaches 2 and applying L'Hopital's Rule. The steady, forced solution 
computed along 45°N is given in Figure 2.3 for a 45° sector of the globe. 
A high pressure ridge is centered over the mountain, which is located at 
the center of the sector, and weak troughs are observed both upstream and 
downstream of the mountain. 

When superposed upon a transient, growing wave, the forced solu- 
tion will cause a cyclonic disturbance to appear more intense both on the 
upstream and downstream side of the mountain, and to appear weaker over 
the mountain. To demonstrate this effect, a developing free Eady-wave of 
approximately the same scale as the forced wave (here, the global wave 
number u, = 8) is taken at various stages of its evolution (Figure 2.4) 
and superposed on the steady solution. The resulting superposition (or 


combined solution) is shown in Figure 2.5. At t = 0, when the free 
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figure 2.3. 


Kee 221.27 Km) 


Steady, forced solution for the linear model presented as 
pg, the deviation from mean-state surface pressure, for a 
450 sector if the globe. Cyclic continuity is assumed at 
the east and west boundaries of the sector. Profile of the 
mountain given by (2.14) is included at bottom. 
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Figure 2.4. Transient, free solution for the linear model given by 
(2.1, 2.5) presented as pg for a 45% sector of the globe. 
Curves are given at six-hourly time intervals at 450N. 


Cyclic continuity is assumed at the east and west boundaries 
of the sector. 
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Figure 2.5. As in Figure 2.4, except for combined (free + forced) linear 
model solution. Profile of the mountain given by (2.14) is 
included at the bottom. 
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and forced solutions are in phase, the combined solution exhibits a 
larger amplitude ridge-trough pattern than the free solution alone. As 
the developing wave moves toward the mountain (at t = 6 and 12 h), the 
amplitude of the combined solution decreases indicating an apparent weak- 
ening, while the amplitude of the free solution continues to grow. As 
the wave reaches the top of the mountain (at t = 18h), it is almost 
completely cancelled in the combined solution by the forced ridge and 
only weak troughs are observed on the upstream and downstream sides of 
the mountain. The cancellation would be exact if the wavelength and 
amplitude of the mountain were the same as, and out of phase with, the 
wave. Between t = 18 and t = 30 h, rapid deepening of the lee-side 
trough is apparent in the combined solution as the free wave moves east 
of the mountain and into phase again with the forced solution. When the 
combined solution is viewed in isolation, it appears as if a well- 
developed cyclonic disturbance dissipates as it approaches the mountain. 
Later, a lee-side disturbance appears to form and amplify rapidly as it 
moves down the mountain slope. 

This behavior is very similar to observations preceding and dur- 
ing Rocky Mountain lee cyclogenesis. In a typical case (Palmen and 
Newton, 1969), a well-developed Pacific cyclone slows and weakens as it 
approaches the west coast of North America. As in the example above, an 
initially strong high pressure ridge centered over the Rockies appears to 
weaken considerably during the next 12-18 hours. However, cyclogenesis 
follows along the lee slopes as rapid deepening and slow eastward move- 
ment are observed during the first 6-12 hours after formation. As the 


cyclone moves away from the lee slopes, the intensification is less rapid 


eZ 


and the eastward phase speed Incas In the words of Palmen and 
Newton, the disturbance "behaves in a manner similar to a cyclone over 
mat terrain." 

A plot (Figure 2.6) of minimum pg versus time for the combined 
and free solutions, shows the major stages of this sequence quite well. 
From t=0 to t= 18 h, the wave in the combined solution appears to 
weaken rapidly (as indicated by rising pg) and suggests a weakening 
surface cyclonic disturbance. Notice that during this period, the free 
solution continues to grow. From t = 18 to t = 36 h, pressure in the 
combined solution falls at a much faster rate than that of the free solu- 
tion, which suggests a more rapid intensification. After t = 36 h, when 
the free wave is away from the mountain, the two solutions appear to 


deepen at the same rate. 





Figure 2.6. Minimum pg versus time at 450N for free (solid) and combined 
(dashed) linear model solutions. 
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The effect of A on the phase of the combined solution 
is shown in Figure 2.7. As the disturbance approaches the mountain from 
t = 6 to 18 h, the phase speed of the combined solution is less than that 
of the free solution. That is, the trough appears to move more slowly as 
observed in Pacific cyclones approaching the windward slopes of the 
Rockies. From t = 24 h, when the lee trough is first observed in the 
combined solution, to t = 36 h, the phase speed of the combined solution 
is less than that of the free solution. This agrees with the observation 
that the lee cyclone moves more slowly during the 6-12 hour period after 
formation over the mountain slope. After t = 36 h, the phase-versus-time 
curves become coincident, which indicate that the disturbance moving away 
from the mountain has a phase speed characteristic of a cyclone over flat 


terrain. 
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Figure 2.7. Phase versus time for free (solid) and combined (dashed) 
linear model solutions. Phase is inferred from the location 
of the minimum value of pe at 450N. 
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It is emphasized that this model is a linear one and does not 
include potentially important nonlinear effects. For example, this model 
requires that the mountain be of a small amplitude so that its effect on 
the mean flow and the developing wave may be ignored. In reality, the 
Rocky Mountain massif does affect the basic airflow. However, the 
success of this conceptual model in describing at least qualitatively the 
major features of lee cyclogenesis indicates the importance of linear 


dynamics in the lee cyclogenesis process. 


C. NON-GEOSTROPHIC MODEL WITH FINITE-AMPLITUDE FORCING 

In the preceding section, the analysis was accomplished with the 
quasi-geostrophic equations and did not consider ageostrophic advections. 
The semi-geostrophic equations described by Hoskins (1975) provide a set 
of equations which is capable of describing nongeostrophic processes such 
as frontogenesis (Hoskins and Bretherton, 1972). This equation set is 
formed by approximating the momentum geostrophically in the equations of 
motion while retaining the complete advecting wind. With a transforma- 
tion of horizontal coordinates, the quasi-geostrophic equations are 
obtained in the new coordinate system. If the quasi-geostrophic equa- 
tions are linear for the problem, these transformed equations are linear 
and may be solved analytically. The solution must then be transformed 
back to physical space for interpretation. The general development which 
follows in section 1 is similar to Bannon (1984) who used the semi- 
geostrophic equations to study the interaction of a front with an 
infinitely long mountain range. In his study, there was no vertical 


Shear in the basic current and the frontogenesis was driven by horizontal 
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deformation. Bannon also used an unbounded atmosphere, while in this 
study, a lid is placed at the non-dimensional height z = 1. 
1. Model Formulation 


As before, consider a Boussinesq fluid. The potential vorticity 


equation may be written in Cartesian coordinates (with z = -ln p/ps) as 
follows: 
dQ =0 ; (2.15) 
dt 
gu ðO dav de dV du, 930 
where = == tat. oe aaa (2.16) 


IZ 3IY 9Z 9X O 3x 3y 92 


The semi-geostrophic approximation is made by replacing (u,v) by (Ug Vg) 
inae GN: 
Using the mean wind of the preceding section, let ug = U = Us + 


Sz, where Us and S are constants. Equation (2.15) may then be re-written 


as 
qq = ‘ (asr) 
dt 
dv dv 2.18 
where A rae ) 
ðZ OX 0 9X ðZ 


Here, the first term on the right-hand side of (2.16) does not appear 
because it is independent of time. If the geostrophic angular momentum, 


M = vg + fx, is introduced, the potential vorticity equation becomes 


36 





IM 39 æM 38  23(M,0 
M , M (M, 0) (2.19) 


MA a 


97 əx  3X əz d(x,z) 


The Hoskins geostrophic coordinate system is introduced with 


v Za 
X= x +2, La Tet ar 
O 
The associated transformation formulas are: 
OR (2.21) 
IX E aX 
A AA (2.22) 
ZAS MONA a 
00 0 
dV f : CES 
where n=f + Í = A ) 
: 1= + 29 
- fọ 9X 
If $, the geopotential in transform space, is defined as 
y? 
b= o> + 9 
2 
where $ is the geopotential in physical space, it follows that 
90 © ð 
DA nO 2a l (act) 
09 ðX 8 ðZ 


O 
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If the initial v is very small, then the potential vorticity will 
be constant, and according to (2.17), it will remain constant. Conse- 


quently, by (2.18), 


39 (2.25) 


The potential vorticity equation may be transformed into the geostrophic 


coordinate system as follows: 





ate 9(M, 0) È 9(M,0) 9(X,Z) oD (M, 0) (2.26) 
9(x,2)  9(X,2) 9(x,2) f Oe) 
and since M = foX, 
30 , (Zn 


= N — 
q 37 


Combining (2.23), (2025). “and (2 2 O lows that 


98 
A 08 foz (2.28) 
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Now, it is convenient to remove the basic stratification by writing 


ò= (z) +Y > (2.30) 
3 gə 
where — = — 
32 0 


Substituting this into (2.29) results in the following 


2 
2y FT 8 gy 
o a f ; MO 


E eee SE 
aX" g To 


which states that the disturbance potential vorticity is zero. The time 


evolution, then, must come from the first law of thermodynamics, 


9 
g- = 0, applied at boundaries z = zs and z - 1. 


The total derivative in the transform space may be written 


0 20 “20 “38 “38 
de ae, (2.32) 
dt a ax y OZ 


Using (2.20), X can be written 


(2.58) 
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where the semi-geostrophic momentum equation has been used to evaluate 
Vg: The first law of thermodynamics now becomes 
30 99 fo - 96 (2.34) 


— + (U. + Sz) — - — Sv +Z—=0 
oT S IX g al 


At the upper boundary, 


A E 


Following Bannon (1984), the lower boundary condition is linearized by 


applying it at Z = 0 rather than z,, and by replacing 


90 90 
37 by the constant 37 


The forced vertical motion is given by 


dh. - hs - hs i (2.36) 


l = — = —— Z 


dt 3X 9X 


Notice that the boundary condition, itself, is not linearized; it becomes 
linear as a result of the assumption that z, = 0. 

Since the interior equation is linear and the boundary conditions 
are linear, the solution may be separated into a steady-state forced 


solution and a free solution: 


Y = Y KoL o (2508) 
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The mountain solution, Yu? satisfies 


2 2 2 
A 
3X? < JA 
g 90/3z 


and the boundary conditions: 


The free solution must satisfy 


2 2 2 
d Ye , me 0 Ye 
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and the boundary conditions: 
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Ss —=0U 


(2.384) 


(2385) 


228 eC } 


(2.39a) 


(22399) 


(22396) 


Notice that the lower boundary condition for the free solution follows 
from setting 7 =Q at Z=0. Both sets of equations are quasi-geostroph- 
ic in transform space. Consequently, the mountain solution is given in 
transform space by (2.13), which gives a high-pressure ridge over the 
crest of the mountain which decays with height. The free solution in 
transform space leads to the Eady (1949) baroclinic instability solutions 
(2.1, 2.5), although the continuous modes (Pedlosky, 1964; Farrell, 1982) 
are required for complete solutions. In transform space, the mountain 
solution and the free solution do not interact. The interaction occurs 
in physical space as expressed through the coordinate transformation 
(20) 
2. Model Application 

Consider the interaction between a small-amplitude, baroclinic 
wave and the flow forced by a finite-amplitude mountain. This requires 
that 

VM VE (2.40) 


and, consequently 
V(X,Z,T) = V(X,Z) = vy (2.4 


The coordinate transformation (2.20) now becomes 


(2.42) 
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The solutions to (2.38) and (2.39) are transformed to physical space with 
(2.42). In this case, the transformation is independent of time because 
it is controlled by the steady mountain solution. 

The effect of the mountain on the amplitude of the disturbance 
may be seen by considering the effect of the transformation on the geo- 
potential. In physical space, the pressure function, $, will change ac- 
cording to 


_ A (2.42) 
eee 


In this application, v = Vy(X,Z), so that (2.43) may be rewritten as: 


Vy (XZ) 
o(X,Z,T) = $(Z) + Oy (X,Z) 3 aa HE a Oe (X,Z,T)-Vy(X,Z)Ve(X,Z,T) 


(2.44) 


This indicates that the primary effect of the onographically forced solu- 
tion is "felt" by the time-dependent, mountain term, Om, because of the 
Squared term. The time-dependent VmMVr term indicates that an orographic 
effect is present. However, there is no orographic effect on amplitude 
at the base of the free trough where Vr = 0. Consequently, changes in 
amplitude (other than those due to the free wave) as the wave moves 
across the mountain are due to the superposition of the forced and free 
solutions. After the wave moves away from the mountain, its amplitude 
will be the same as if there had been no mountain. In other words, for 
the assumptions made above, a finite-amplitude mountain has no effect on 


the growth of a small-amplitude baroclinic wave. Even over the 
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mountain, the minimum ee (where Pe is the time-dependent part of Ê) is 
the same whether or not it is over the mountain. Only its location will 
be affected by the mountain. 

The effect of the mountain on the phase of the disturbance may be 


seen by diffentiating (2.42) with respect to X: 


ð X 2 ee (2.45) 
IX f ox 


This result indicates that the horizontal dimension of the solution will 


be stretched in physical space when 


ay Jy 
ae O and compressed when ay 0 


In physical terms, the stretching of the horizontal dimension corresponds 
to a phase acceleration of the solution, and the compression represents a 
deceleration. These effects are caused by Uy, the divergent component of 
the wind, which is forced by the mountain. The divergent component Un, 
iS a maximum at the mountain top because there is divergence on the up- 
wind side and convergence on the downwind side. Consider Figure 2.8 in 
which the solution Vm(X,0) is plotted in transform space for a mountain 
profile H(X). The mountain forces a northward deflection of the wind 
(Vy 7 0) on the windward side and a southward deflection on the leeward 
Side. If -X, and +X, represent the points at which 


dy 
a 0, a wave will be decelerated when it is west of -X1, 
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Figure 2.8. Steady, forced solution of the semi-geostrophic poded , VM: 
in transform space. At the points -X1 and +X], “Wm/"X = 0; 
these points separate the region of phase acceleration from 
the regions of deceleration. 


accelerated when it is between -Xi and +X¡ (over the mountain), and de- 
celerated again when it is east of +X}. 

The phase acceleration may be seen by comparing the free 
solution (of (2.38)) in transform space and physical space. As indicated 
previously, the free solution is given by (2.1, 2.5). It is transformed 
to physical space by reversing the Hoskins transformation to obtain 


LT i T 
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where Vm is computed by taking the derivative of (2.9, 2.13) with respect 
to X. The solutions in transform space and physical space may then be 
computed using the same constants as in section II.B.2, except 
Zags 3.0 km. The result is shown in Figure 2.9. The horizontal scale of 
the developing Eady wave is approximately the same as the scale of the 
mountain (global wave number W, = 8). To- focus on the acceleration 
caused by the mountain, the trough in transform space is initially posi- 
tioned in phase with the trough in physical space at the point X = -X1. 
This point where the acceleration effect begins is determined analytical- 
ly using (2.42). In this way, the phase difference between the two 
troughs on the lee side of the mountain will be a maximum. 

As the wave begins to move up the mountsin slope (t = 6 h), evi- 
dence of acceleration is already present in the physical solution. At 
t = 12 h, when the trough in transform space reaches the top of the 
mountain, the trough in physical space is located half-way down the lee 
Side of the mountain. At t = 18h, the trough in transform space has 
moved over the lee slope and now trails the trough in physical space by 
nearly 500 km. This phase difference has decreased by t = 24 h, which 
indicates that the trough in physical space has slowed. This effect on 
phase speed is caused by the divergent component of the wind which is 
forced by the mountain. Bannon (1984) obtained similar results with a 
front moving over a mountain range. After the disturbance moves away 
from the mountain and Vm(X,Z)”0, the phase speed of the physical solution 
1s approximately the same as that of the transform solution, and the lee 


disturbance now behaves similar to a cyclone over flat terrain. 
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Figure 2.9. Semi-geostrophic model solutions presented as p< along 450N 
in transform space (solid) and in physical space (dashed). 
The 3.0km mountain is included at the bottom. 
This latter result is similar to that discussed in the preceding 
section, but, in this case, the acceleration and deceleration of the dis- 


turbance is actually caused by the interaction with the mountain. As the 


disturbance moves down the lee slope and away from the mountain it is 


47 


decelerated and moves more slowly than a similar disturbance over flat 
terrain. After the disturbance has left the mountain's influence, its 
phase speed increases to that of the disturbance over flat terrain. As 
in the preceding section, a relatively simple model has been used to de- 
scribe key aspects which are observed in Rocky Mountain lee cyclogenesis. 
The results of this chapter suggest that the cyclogenesis which 
is observed in the lee of the Rocky Mountains is largely unaffected 
dynamically by the presence of the mountain. In both models, the rapid 
intensification which is characteristic of lee cyclogenesis is observed 
as the disturbance moves down the lee slope of the mountain. However, in 
neither model is there a net orographic effect on the amplitude of the 
disturbance. The characteristic phase da of cyclonic disturbances 
on both the windward and leeward sides of the mountain is also observed 
in these models. In the linear model, the effect was due to superposi- 
tion of two solutions; in the semi-geostrophic model, the larger mountain 
affects the phase speed of the wave. However, after the disturbance has 
moved away from the mountain, it shows no net orographic effect. The 
success of these models in explaining the major features of Rocky Moun- 
tain lee cyclogenesis suggest that it is a manifestation of simple bara- 
clinic instability distorted by the mountain - and not any new or 
enhanced instability mechanism. This idea will be tested in Chapters V 
and VI where the primitive equations are integrated with more realistic 


initial conditions. 
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III. ATMOSPHERIC MODELS 


Two atmospheric prediction models are used in this study: the UCLA 
general circulation model described by Arakawa and Lamb (1981) and the 
Naval Environmental Prediction Research Facility (NEPRF) global spectral 
model developed by Dr. T. Rosmond. These two models represent state-of- 
the-art atmospheric research models Fa have many desirable features. 

The UCLA model is designed specifically to provide an accurate repre- 
sentation of air flow over and near topography. The basic numerical 
scheme was developed for the shallow water equations which conserve abso- 
lute potential vorticity. The UCLA finite-difference scheme is designed 
to conserve the domain-averaged potential vorticity while it approximate- 
ly conserves the potential vorticity of individual parcels. Arakawa and 
Lamb (1981) have shown that, for a given coarse grid, simulations of the 
airflow over steep topography improved significantly when this scheme is 
used. In this study, airflow over a long, narrow mountain is considered, 
and it is believed that the UCLA model provides the best available 
representation. 

The NEPRF model es a spectral version of the 1977 UCLA model. Al- 
though it is not explicitly designed to conserve the domain-averaged 
potential vorticity, it is quite accurate because of the spectral repre- 
sentation. It is used in this study to investigate the continuous-mode 
growth mechanism. Because this mechanism is a linear one, a linearized 
model is sufficient for the investigation. The NEPRF model was selected 


primarily because it can be linearized more easily and run more 
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rapidly than a comparable finite-difference model. Since the physics are 
very similar to those of the UCLA model, the results should be about the 
same. Additionally, the fully nonlinear version of the NEPRF model can 
be used as a check on the effect of truncation errors in the finite- 
difference runs. The following sections describe these two models and 


the integration schemes used. 


A. THE UCLA MODEL 

The UCLA general circulation model is described in detail by Arakawa 
and Lamb (1977, 1981). A short summary of the version used in this study 
is provided here. The model consists of the primitive equations for an 
inviscid, adiabatic, and hydrostatic atmosphere; moisture and its effects 
are not included. The prognostic variables are the horizontal components 
(u, v) of the wind velocity, potential temperature (8), and pressure 


(7 = po - pt). The model's vertical coordinate is O, which is defined as 


S , (3.1) 
Poga Pt 


S 
where ps iS surface pressure and pt is the top of the model atmosphere. 
In this study, pt = 200 mb. Channey and Drazin (1961) have shown that 
wavelengths shorter than planetary scale are trapped in the vertical. In 
addition, the results of Section II.B (Figure 2.2) indicate that the ono- 
graphic wave is strongly damped in the vertical. Consequently, model 


solutions should not be significantly affected by the boundary at 200 mb. 
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The domain is a 450 sector of the Northern Hemisphere with a wall at 
the equator and cyclic continuity at the east and west boundaries. The 
model atmosphere extends from the earth's surface to 200 mb, and 4 ¡s 
assumed to vanish at both of these boundaries. The variables are stag- 
gered horizontally according to Arakawa's Scheme C (Figure 3.1) with a 
grid spacing of approximately 2.80 longitude by 2.750 latitude. Vari- 
ables are staggered vertically as shown in Figure 3.2 in six layers 
spaced equally in%. Rennick and Williams (1985) have shown that with an 
upper boundary at 200 mb, the equal spacing in % gives well-behaved 
results for the forced mountain solution at the horizontal scales con- 
sidered in this study. Spatial derivatives are approximated using a 
E inorder finite difference scheme. The vertical differencing scheme, 
which was developed by Arakawa and Suarez (1983), has excellent integral 
properties and it eliminates the systematic error in the hydrostatic 
equation which was present in previous models. A nonlinear horizontal 
diffusion of momentum is included to parameterize the effects of subgrid- 


scale motions. 
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Figure 3.1. Horizontal distribution of variables in the UCLA model grid. 
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Figure 3.2. Vertical distribution of variables in the six-layer version 
of the sigma-coordinate system. 


The model integration proceeds in a series of Euler-backward time 
steps followed by five centered time steps. Convergence of the meridians 
toward the poles would normally require the use of an extremely short 
time step to maintain computational stability. To avoid this require- 
ment, the zonal smoothing technique described by Arakawa and Lamb (1977) 
is used. In this technique a local stability criterion, S, is used; it 


is defined as follows: 


S = (AA/d>)cos 6/ sin(ndA/2) (At, /At) (Sua 


where AA = longitudinal grid size, 4d = latitudinal grid size, 


ât = 360 s, At = time step, and n = wave number. 
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Smoothing is performed only when S < 1; this occurs generally at high 
latitudes and for higher wave numbers. The smoothing is accomplished by 
reducing the amplitudes of the Fourier components of the zonal pressure 
gradient and the zonal mass flux by the factor S. As indicated by (3.2), 
the amount of smoothing may be reduced by decreasing the time step. To 
measure the effect of the zonal smoothing on the features of interest in 
this study, runs were made in which a westerly mean flow was forced over 
a 1.5 km mountain using three-, six-, and ten-minute time steps. Compar- 
ison of the 36-h forecasts of surface pressure and 500 mb vorticity 
revealed virtually no differences between the solutions with three- and 
six-minute time steps. Therefore, it was concluded that a six-minute 


time step was adequate. 


B. THE NEPRF SPECTRAL MODEL 

A linearized form of the NEPRF spectral model is used in this study 
to investigate the continuous-mode growth mechanism. A short summary of 
the basic model and the version used in this study is presented in this 
section. 

The spectral formulation was accomplished by Dr. T. Rosmond of NEPRF 
and is described by Lubeck, Rosmond and Williams (1977). The model's 
basic equations are obtained from the nonlinear primitive equations for 
an inviscid, adiabatic, and hydrostatic atmosphere. the prognostic vari- 
ables are the horizontal components of the wind velocity (u, v); temper- 
ature (T); and the natural log of terrain pressure (q = In pc). The 
variables are represented spectrally in the horizontal with wave number 
40 triangular truncation. The model's vertical coordinate is 9, which is 


given by (3.1), with py = 0. The vertical derivatives are approximated 
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using finite differences. Vertical velocity, o, is assumed to vanish at 
the upper and lower boundaries. The model employs a semi-implicit time- 
differencing scheme which allows a 15-minute time step. 

The version of the model used in this study is linearized according 
to Rennick and Williams (1985). This is accomplished by forcing the time 
Paden to be zero at all wave numbers except the one of interest. 
The domain is restricted to one wavelength in the east-west direction, 
and cyclic continuity is assumed. To simulate the wall placed at the 
equator in the UCLA model, the Northern Hemisphere conditions are mirror- 
ed in the Southern Hemisphere. Six- and 20-layer versions of the model 
are used. A comparison of the results of these two versions will give an 
indication of the effect of the vertical resolution. The layers are 
spaced equally in o. Rennick and Williams (1985) found that this spacing 
yields well-behaved solutions. Because the experiments are limited to 
the linear domain, the use of semi-implicit time-differencing allows a 


30-minute time step. 
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IV. INITIAL CONDITIONS AND TOPOGRAPHY 


Many factors must be considered in the development of the initial 
conditions for experiments which treat the lee cyclogenesis problem. The 
wind shear (both vertical and horizontal) of the mean flow determines its 
Stability and, consequently, plays a controlling role in the development 
of the cyclone. In addition, the vertical structure of the disturbance 
is critical to the rate at which mean kinetic energy is converted into 
eddy kinetic energy in a developing wave. Finally, as indicated by the 
climatology of cyclogenesis frequency, the size, scale and orientation of 
topography appear to play a role in the occurrence of cyclogenesis. Con- 
sequently, the initial conditions must be selected with care. 

In this study, an attempt is made to reproduce the observed charac- 
teristics of lee cyclogenesis as closely as possible using analytical 
expressions. The structure of the mean flow is similar to that presented 
in Palmen and Newton (1969) in their description of cyclogenesis. The 
disturbance is specified using the results of Chapter II. The topography 
is patterned roughly after the Rocky Mountains. In the following sec- 
tions, the specific analytical expressions and techniques used to initia- 


lize the model are described. 


A. MEAN FLOW 

In this study, lee cyclogenesis is viewed as a mid-latitude manifes- 
tation of baroclinic instability. Consequently, the initial mean state 
for all experiments consists of a baroclinically unstable westerly 


current. The specific form of the upper-level wind is 
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= - 2 > (4.1a) 
Uyló ,p(S,9 )] = Upsech [Y (9 - $ 3] (1n(ps/p)/1In(ps/Pmax) 


0 i (4.1b) 


TEE 


where Ug = 40 m So aor = 459N, po = 1013.25 mb, pmax = 200 mb, and Y is 
the halfwidth of the jet. The maximum wind at each level occurs at 45°N. 
Wind speed varies linearly with the log of pressure and approximates 
profiles given by Palmen and Newton (1969) in case studies of 
cyclogenesis. The latitudinal variation is taken from Haltiner and 
Williams (1980). By setting the halfwidth of the jet, Y, equal to 89 
latitude, a horizontal profile which agrees quite well with that observed 
prior to cyclogenesis results. However, this profile meets the necessary 
condition for barotropic instability in regions to the north and south of 
the jet stream. Consequently, a profile in which Y = 169, that does not 
meet this criterion is included in this study to isolate the baroclinic 
effects. Meridional cross-sections of the two initial mean state 
velocity fields are shown in Figure 4.1. 

A mean surface current is included with the same latitudinal struc- 


ture as the upper-level flow, 


uc(%) = Uggsech LY(> - >5)] > (4.2a) 


(4.2b) 
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Figure 4.1. Meridional cross-sections of the initial mean wind, u, for: 
a) 169 jet, and b) 89 jet. Contour interval is 10 m s”?. 


where y and @ are as defined above and Ugg = 5m s-+. The upper-level 
wind given by (4.1) is modified by adding the surface current and apply- 


ing the gradient correction to account for the earth's sphericity: 


u[o,p(o,6)] = -Ra cos > + 


fa cos 6[1+(2(u + u )/na cos 9)]0-5 
u S 


where a is the radius of the earth, 2 is the earth's rotation rate, and 


Uy and us are given by (4.la) and (4.2a), respectively. 


oy 


Mean surface pressure and surface temperature are initialized from 


the mean surface wind using the geostrophic wind equation: 


z p $ : 
ps($) = expLIn pc(d9) - (a/R) J fus(9)/Ts(9) de] > (4.4) 


O 
where Ps(þo = 450N) = 1013.25 mb, us(ẹ) is given by (4.2a), and Ts (¢) 
represents the mean surface temperature. The integral in (4.4) is evalu- 
ated using a Simpson's rule approximation. Because the latitudinal 
structure of To is unknown, ps(?) is computed interatively as follows: 
an initial guess is made for Te: ps(9) is computed using (4.4); Te is 
then adjusted to the NACA standard temperature for the computed De and a 
new Ps($) is derived; this latter step is repeated until the adjustment 
of Te is less than 0.01°K. The solution converges in approximately 10 
iterations. 

Upper-level temperature is specified by integrating the geostrophic 


thermal wind equation which gives: 


: : : (4.5) 
T(o,p) = T(o5,p) + (a/R) J f(3u,,/Tn p) de 


O 


where Ts(do,P) is the temperature at 45%N and u, is given by (4.1). The 


integral in (4.5) is also evaluated using Simpson's rule. 
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These initial conditions were integrated over a flat lower boundary 
to test the balance of the mean state. Harmonic analysis of the v-com- 
ponent of the wind at all levels and at a variety of latitudes showed 


wave number zero amplitude fluctuations of less than 1 m oe 


B. DISTURBANCE 

Two types of disturbances are used in this study: a barotropic dis- 
turbance (amplitude and phase independent of height), and an orographic- 
ally-forced disturbance which decays exponentially with height. The 
barotropic disturbance is used in the majority of the experiments to 
Study the effect of topography on the evolution of a baroclinic wave. 
The orographically-forced disturbance is used only in Chapter V to study 
the continuous-mode growth mechanism. 

The barotropic disturbance consists of a weak wave which varies sinu- 
soidally with longitude. The maximum amplitude of the disturbance occurs 
at 45°N. Fields are balanced geostrophically with a constaft f and are 


given by 


®' = fo A*sin(nA)sin*(20) i (476a) 
p' = p®'/RT A (4.6b) 
u' = -(1/fya) 99'/9b ; (4.6c) 
v' = 1/fya cos 0) 909'/9A s (4.6d) 
O à (4.6e) 


where T = 2730K and p = 1013.25 mb. Test integrations were conducted in 


which the initial wave number of the disturbance was varied to determine 


BY 


the most unstable wavelength. The wave number 8 perturbation exhibited 
the most rapid growth and is the one used in this study. 

The orographically-forced disturbance is specified using the steady 
forced solution to Eady's model given by (2.9) and 2.13). The fields are 


balanced geostrophically: 


Ie i ; (4.7a) 
pg = -p/RT[(Wshse 2 COS yx sin2(24))/(H(Usog  - Sete )] 
ay Ushs c% Z 
US is di o psinh -r + 
af HU or=Se YE? £ 
2 4.7b 
yecosh Ecos ux (2 sin 49) \ 
= 
] Ushs e+y z 
v! = 2. = | S mit 
9 HS € 
E HU or Se YE? € 
Z : , (4.70) 
ypcosh 4 | cos uX (2 sin 22 y) 
1 99' TUshs Z 
Ths Se» > | cosh- t 
3 RH(U. ofS e YE) e” 
a (4.7d) 


YE sim E cos ¡Xx (sin? g) 


where — = 1.167 x 10%més-2, Us is the mean surface wind, hs is the ter- 


rain elevation, H = 8 km is the scale height of the atmosphere, and z is 


non-dimensional; height (here Z= 1 - o). 
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C. TERRAIN 
The topography used in this study is designed to resemble the Rocky 


Mountains as a long, meridional barrier to westerly flow. It is given by 





A-A 
ZO EOS I =) a » JA =- r < 4dr (4.8a) 


0 A AN 


where AA is the longitudinal grid spacing and Ag is the longitude at 


which the mountain is centered; Z (9) is given by 





Zo EE A 
o-> , 
Le cos? pe ‘| 3 On + 3A > > > On 
344 2 
2 (0) = (4.8b) 
LS EA ets, i » 5 > > > ds - 349 
3Ao 2 


0 , elsewhere 


where Zo is the mountain height, Ad is the latitudinal grid spacing, ón 
is approximately 61.75%N and bs is approximately 31.250N. The resulting 
mountain (Figure 4.2) is 22.50 wide and extends from approximately 230N 


to 700N. 
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Figure 4.2. Plan-view of the topography on a 45% sector of the globe. 
Vertical line at the center of the figure represents the 
mountain ridge line. 

Terrain was initialized in the model by raising the mountain from 
zero to the values in (4.8) during the first 12 hours of the integration. 
This technique is designed to minimize the generation of inertia-gravity 
waves during the initial part of the integration while the mean flow, 
which is analytically balanced without topography, adjusts to the pre- 
sence of the mountain. Topography is incremented at each time step 


according to: 


Tt 
2 (1,9) sinf (=) „t< 12 hours ie 
ik he | m 24 (4.9) 
m 


z (A,¢) , t > 12 hours 
m* 


This technique has been used successfully by Tibaldi et al. (1980), 
Williams et al. (1981), and Walker (1982) and has been found to produce a 


smooth and rapid adjustment (within approximately 6-12 hours). 


62 


V. CONTINUOUS-MODE GROWTH EXPERIMENTS 


The results of Chapter II suggest that lee cyclogenesis may be inter- 
preted as a manifestation of ordinary baroclinic instability distorted by 
the mountain. In this chapter, the suggestion by Farrell (1982) that 
continuous-mode growth plays a significant role in lee cyclogenesis is 
considered. As indicated in Chapter I, the initial conditions are 
crucial to this mechanism and must be based on a realistic flow. In this 
Study, the mean flow is assumed to contain no perturbations initially, 
and, at the same time, t, a disturbance composed of a spectrum of dis- 
crete and continuous modes is excited in response to orographic forcing. 
This might be expected if the cross-topography flow increased rapidly 
from zero. Two sets of numerical experiments are conducted in which an 
orographically-forced disturbance is allowed to evolve for at least 24 
hours in a baroclinically unstable mean flow. In the first set, a linear 
model is used to examine the evolution of an orographically-forced dis- 
turbance which is specified initially using the forced Eady solution 
described in section II.B. Use of the linear model is justified because 
the continuous-mode growth described by Farrell is a linear mechanism. 
In the second set of experiments, the possibility is considered that the 
Structure of the orographically-forced disturbance is more complex than 
the forced Eady solution. In these experiments, the forced disturbance 
is not included in the initial conditions but develops during the inte- 


gration in response to the mean flow being forced over the topography. 
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In addition, a nonlinear numerical model is used to include all possible 
orographic effects on the evolution of the continuous modes. 

Results of these numerical experiments are examined for evidence of 
enhanced growth which could indicate a significant contribution to the 
amplification of the disturbance. According to Farrell's analysis, con- 
tinuous-mode growth could occur for approximately 6-8 hours after excita- 
tion for the magnitude of vertical shear used in these experiments. 
Pedlosky (1964) has shown that, in the long term, these modes damp as 
1/t2. Consequently, if the rapid algebraic growth predicted by Farrell 
is present, it should be evident within the first 24 hours of the 


integration. 


A. EVOLUTION OF VERTICALLY TRAPPED DISTURBANCES IN A LINEAR MODEL 

In this section, the excitation of an orographically-forced, vertic- 
ally-trapped disturbance is considered as a possible trigger mechanism 
for rapid continuous mode growth. Numerical experiments are carried out 
in which the forced Eady solution described in section II.B is allowed to 
evolve in a baroclinically-unstable, zonal mean flow. Topography is not 
included because the cont inuous-mode-growth mechanism is linear, and, 
Once excited, these modes grow independently of the mountain. In any 
case, the mountain solution is independent of time and it cannot interact 
with the disturbance since the model equations have been linearized. 

The six-layer, linearized version of the NEPRF spectral model is used 
to carry out the experiments. The mean state is given by (4.3-4.5) with 
y = 8°, and the initial disturbance is given by (4.7a-d). The distur- 


bance is allowed to evolve for at least 24 hours. Model runs are 
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carried out with global wave number LH, = 8, 12, and 16 because these wave 
numbers represent disturbance wavelengths of interest in studies of lee 
cyclogenesis. 

The results of the wave number 8 experiment are shown in Figure 5.1. 
After an intial period in which the disturbance amplitude damps, amplifi- 
cation is observed at © = 0.9167 (p = 930 mb) between t = 6 and 12 h. 
However, this growth is confined to the lowest level of the model and is 
not sustained beyong t = 12 hours. Little if any growth is observed at 
upper levels. The model run is extended to t = 60 h to see if the dis- 
crete mode growth predicted by baroclinic instability theory is observed. 
In Figure 5.2, the amplitude of v at U, = 8 for level 6 is plotted versus 
time. Except for the growth spurt between t = 6 and 12 h noted above, 
Virtually no growth is observed during the first 42 hours of the run. 
However, beginning at approximately t = 45 h, growth of the most unstable 
discrete mode is indicated by the rapid amplification. 

The results of the wave number 12 and 16 experiments are shown in 
Figures 5.3 and 5.4, respectively. Behavior in both cases is roughly 
Similar to that observed in the wave number 8 experiment. A small spurt 
of growth at the lowest level is observed in Figure 5.3 between t = 12 
and 18 h between t = 6 and 12h in Figure 5.4. Again, little if any 
growth occurs at higher levels. 

Walker (1982) suggested that higher vertical resolution may be 
necessary to model the dynamics of lee cyclogenesis properly. Farrell 
(1982) has also shown that continuous-mode growth is favored in distur- 
bances with higher wave numbers, which indicates that a higher vertical 


resolution may be necessary for the model. The three experiments 
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Amplitude of wave number 8 of v at 45%N as a function of 9 
for linear Experiment 1. Initial amplitude (dashed) versus 


o is included. 
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Amplitude of wave number 8 of v at o = 0.9167 at 450N 


Pigúure 5.2. 
plotted as a function of time for linear Experiment 1. 
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Figure 5.3. Amplitude of wave number 8 of v at 45°N as a function of o 
for linear Experiment 2. Initial amplitude (dashed) versus 
o is included. 





Figure 5.4. Amplitude of wave number 8 of v at 450N as a function of o 


for linear Experiment 3. Initial amplitude (dashed) versus 
o is included. 
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described above were re-run using a 20-layer version of the NEPRF model 
to determine if more than six layers is required to resolve continuous- 
mode processes. In Figures 5.5a, b, and c, cross-sections at 459N are 
shown for the 20-layer solutions superposed on the six-layer solutions at 
tale fh. In Figure 9.5a, the solutions are nearly coincident at all 
levels, which indicates that there are essentially no differences between 
the six-layer and 20-layer model runs. In Figures 5.5b and c, anp AERE 
differences between the two solutions are only slightly greater (less 
than lm sy and the phases are very close. No evidence of rapid 
growth was observed in any of the high resolution runs. Therefore, it is 


concluded that the six-layer resolution is sufficient for this research. 


B. DISTURBANCED FORCED BY THE TOPOGRAPHY IN A NONLINEAR MODEL 

In the preceding section, the structure of the forced wave was expli- 
citly specified in the initial conditions and a linearized model was 
used. The possibility exists that a more complicated structure is in- 
volved than given in the analytically-specified initial condition. In 
this section, cases are examined in which no attempt was made to control 
the initial structure of the forced orographic wave. Numerical experi- 
ments are conducted using the six-layer UCLA model, described in Chapter 
III, which includes nonlinear processes. The initial conditions consist- 
ed solely of the baroclinically unstable mean flow given by (4.3-4.5). 
In each experiment, the mean flow is forced over a long, narrow mountain 
(given by (4.8)) and was allowed to evolve for 90 hours. 

Four experiments were conducted. In the first, the half-width of the 
jet stream, Y, was 80 and the mountain height, Zm,» was 1.5 km. In the 


second experiment, Y = 80 and Zm, = 3-0 km were used. In the third 
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experiment, Y = 160 and Zm, = 1.5 km, and in the fourth, Y = 169 and 
Zm, = 3.0 km. Topography was raised from zero to its final height in 
12 hours according to (4.9). Forcing of the mean flow by the mountain 
theoretically will excite a spectrum of discrete and continuous modes. 
As in the preceding section, for continuous modes to be significant, 
growth must occur earlier and be more rapid than that predicted for the 
discrete modes. 

Results of the first experiment (Y = 80, Zm, = 1.5 km) are shown in 
Figure 5.6. When the mountain reaches its final height at t = 12 h, the 
forces ridge-trough pattern predicted by the Eady solution is observed. 
However, the high-pressure center is located over the lee A of the 
mountain rather than over the ridge line as in the results of section 
II.B. The high-pressure ridge is maintained over the southern portion of 
the mountain throughout the model integration. However, in the region 
beneath the jet stream, the eastward elongation of the ridge indicates 
that a portion of the high is advected downstream. These results are in 
qualitative agreement with those of Huppert and Bryan (1975) who found 
that topographically forced circulations were advected downstream in 
strong mean currents and were bound to the topography in weak mean flows. 

Examination of this sequence does not reveal any rapid continuous- 
mode growth. Between t = 12 and t = 24 h, when the continuous-mode 
growth should be evident, there is no significant cyclonic development. 
A weak trough appears at t = 24 h, but only slight amplification is 
observed during the following 24-30 hour period. This is certainly not 


the type of growth normally associated with lee cyclogenesis. 
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Figure 5.6a. Sea-level pressure contours for nonlinear Experiment 1 at 
12 h. Contour interval is 2.5 mb. Domain is a 45° sector 
of the globe with cyclic continuity at the east and west 
boundaries. 





l 


Figure 5.6b. Sea-level pressure contours for nonlinear Experiment 1 at 
18 h. Contour interval is 2.5 mb. Domain is a 450 sector 
of the globe with cyclic continuity at the east and west 
boundaries. 
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Figure 5.6c. 


Figure 5.6d. 





Sea-level pressure contours for nonlinear Experiment 1 at 
24 h. Contour interval is 2.5 mb. Domain is a 45% sector 
of the globe with cyclic continuity at the east and west 
boundaries. 





Sea-level pressure contours for nonlinear Experiment 1 at 
30 h. Contour interval is 2.5 mb. Domain is a 45° sector 
of the globe with cyclic continuity at the east and west 
boundaries. 


72 


Figure 5.6e. 


Figure 5.6f. 





Sea-leve] pressure contours for nonlinear Experiment 1 at 
36 h. Contour interval is 2.5 mb. Domain is a 450 sector 
of the globe with cyclic continuity at the east and west 
boundaries. 





Sea-level pressure contours for nonlinear Experiment 1 at 
42 h. Contour interval is 2.5 mb. Domain is a 45% sector 
of the globe with cyclic continuity at the east and west 
boundaries. 
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Figure 5.6g. Sea-level pressure contours for nonlinear Experiment 1 at 
48 h. Contour interval is 2.5 mb. Domain is a 45% sector 
of the globe with cyclic continuity at the east and west 
boundaries. 


In Figure 5.7, the amplitude of v for global wave number 8 at 
o = 0.9167 is plotted versus time. In the first 12 hours, the amplitude 
grows rapidly in response to the "growing" mountain. Between t = 12 and 
t = 24 h, the model atmosphere adjusts to the presence of the mountain. 
A weak inertia-gravity wave (amplitude approximately 1-2 m s71) is ob- 
served. However, this wave does not grow and, in fact, appears to damp 
slightly with time. Amplification is not observed between t = 12 and 
t = 50 h. After t = 50 h, wave number 8 begins to amplify at an approxi- 
mately exponential rate. This type of growth is predicted for distur- 
bances in baroclinically unstable mean flows, but it is normally associ- 


ated with the most unstable discrete mode, and not with continuous modes. 
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Figure 5.7. Amplitude of wave number 8 at o = 0.9167 at 450N versus time 
for nonlinear Experiments 1 (solid) and 2 (dashed). 


Results of the second experiment (y = 89, Zm, = 3.0 km) are similar 
as seen in Figure 5.7. As expected, the larger mountain forces a 
Stronger deflection of the mean flow during the first 12 hours. Doubling 
of the mountain height approximately doubles the amplitude of the re- 
sponse. A more intense forced ridge-trough pattern is observed in the 
sea-level pressure contours (not shown), although no evidence of rapid 
growth is observed during the first 24 hours of the integration. Growth 
of the most unstable discrete mode is indicated here also, and it begins 


at about the same time (t = 50 h). 
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Results (Figure 5.8) of the third (y = 169, z US km) and fourth 
(y = 160, Zm, = 3.0 km) experiments are qualitatively similar. The 
response to the growing mountain during the first 12 hours is approxi- 
mately the same ad in the y = 89 runs, and doubling the mountain height 
approximately doubles the amplitude. A smaller gravity-wave oscillation 
(amplitude less than 1 m a is observed after t = 12 h, which indicates 
a better geostrophic adjustment. No initial growth which might indicate 
growing continuous modes is observed. In these two runs, baroclinic 
growth associated with the most unstable discrete mode is observed after 


t = 40 h (about 10 hours earlier than in the y = 80 runs). 





Figure 5.8. Amplitude of wave number 8 at o = 0.9167 at 450N versus time 
for nonlinear Experiments 3 (solid) and 4 (dashed). 
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The results of this chapter do not support the hypothesis that con- 
tinuous-mode growth plays a significant role in lee cyclogenesis. Evi- 
dence of sustained, rapid growth is not observed during the first 24 h of 
any of the numerical experiments examined in this chapter. In the linear 
experiments, the vertically-trapped disturbance shows spurts of growth in 
the lower levels, but this growth is not sustained beyone approximately 
3 h after onset. The overall character of these solutions appears to be 
oscillatory rather than amplifying. | 

In the nonlinear experiments, slight amplification of the forced sur- 
face-pressure trough is observed after 12 hours. However, since no evi- 
dence of rapid initial growth is observed in the v-component of the wind 
at 0 = 0.9167, amplification must be confined to a region close to the 
earth's surface. The growth observed in these experiments is attributed 
to the most unstable discrete mode - and not to continuous modes. Con- 
sequently, if as Farrell suggests, growing continuous modes play a role 


in lee cyclogenesis, other justifiable initial conditions must be found. 
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VI. BASIC NUMERICAL EXPERIMENTS 


In this chapter, enhanced baroclinic instability and superposition 
are considered as possible lee cyclogenesis mechanisms. Markine's (1975) 
analysis of quasi-geostrophic flow past an infinitely long mountain 
Showed increased vertical wind shear on the lee side of the mountain, 
which indicates enhanced baroclinicity. He has suggested that as a dis- 
turbance moves into the lee of a mountain complex, it encounters a 
Signficantly destabilized mean flow and rapid cyclogenesis results. 

Merkine's suggestion is investigated in this chapter using a series 
of numerical experiments in which a free disturbance evolves in a baro- 
Clinically unstable flow. To isolate the effects of the mountain, runs 
are made with and without topography. One of the two baroclinically un- 
stable mean states that are considered also meets the necessary condition 
for barotropic instability. Mountain height and length, as well as the 
initial amplitude of the disturbance, are varied. 

Results are examined for evidence of increased growth which could 
indicate enhanced instability on the lee side of the mountain. As indi- 
cated in Chapter Il, observed lee-side growth could be explained by a 
superposition of the growing wave on a steady, orograpnically forced lee- 
side trough. Consequently, for the enhanced-instability mechanism to be 
considered a significant contributer to lee cyclogenesis, it must be 
demonstrated that disturbances moving over mountains develop more rapidly 


than disturbances moving over flat terrain. That is, it must be shown 
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that this development significantly exceeds that caused by a combination 


of superposition and normal baroclinic growth. 


A. NUMERICAL EXPERIMENTS 

Numerical experiments are carried out using the version of the UCLA 
model described in Chapter III. The initial mean state consists of a 
baroclinically unstable, westerly current balanced analytically for a 
flat lower boundary (that is, no topography). The basic current, given 
by (4.1-4.5), possesses both horizontal and vertical shear and reaches a 
maximum speed of 45 m s-1 in the top layer of the model at 450N. A mean 
westerly surface current with a similar horizontal profile as the upper 
level flow is included. The maximum surface wind is 5 ms7! at 450N. 
The mountain consists of a long, smooth barrier (approximately 47° lati- 
tude by 249 longitude) oriented north-south as described by (4.8). 

1. Experimental Techniques | 

Each experiment consists of three runs: Control, Interaction, and 

Non-interaction. In the Control run, a weak disturbance, given initially 
by (4.6), is allowed to evolve into a mature cyclone with no topography 
included. The time-dependent solution for this run may be represented 


as: 
Welt) = Yp + Ye (t) 
where Ye is a zonally symmetric, time-independent state in gradient 


balance over flat topography. Here, A is the instantaneous deviation of 


the solution from the mean state. Because the UCLA model is fully 
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nonlinear, there are no implied constraints on the magnitude or mean 
value of ¥ A(t). 

In the Interaction run, a growing baroclinic wave evolves as it 
passes over the mountain. The initial conditions for the integration are 
constructed by combining r (taken from the Control run at some arbrit- 
rary time ty using (6.1) with an orographically forced, mountain mean 


state such that 


Y(to) = YM + Y¿(to) 


Fm is a zonally asymmetric, time-independent mean state representing flow 
in the presence of the mountain. These conditions are integrated forward 
in time with full-height topography until the disturbance passes over the 
mountain. The resulting Interaction solution may be formally separated 


into two parts analogously to (6.1): 


Y] (t) = Yy + Y5(t) 


To isolate the degree of interaction between the growing wave and the 
mountain-induced distortions in the mean state, a Non-interaction solu- 
tion ¥y(t) is constructed without model integration by combining Yalt), 
taken from the Control run at each desired output time, with ym such that 


- : (6.4) 


wy(t) = WM t we(t) 
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In this .way, the Non-interaction solution parallels the development in 
Chapter II in which the time-dependent Eady solution was superposed on 
the steady, orographically forced solution to describe the passage of a 
wave over the mountain. 

The observed growth is affected by the superposition of the 
disturbance on the forced lee-side trough in both the Interaction and 
Non-interaction runs. However, the growth of the disturbance in the 
Interaction run is also affected by the stability of the orographically 
forced mean flow, whereas the growth of the Non-interaction disturbance 
is affected by the stability of a zonally uniform mean flow. Therefore, 
comparison of these two runs will reveal differences in growth caused by 
the orographically forced distortions to the mean flow. If mean flow 
destabilization is a significant factor in lee cyclogenesis, large 
increases in growth should be observed in the Interaction runs. 

2. Derivation of the Mountain Mean State 

The key to the experimental procedure described in the preceding 
section is the steadiness of the mountain mean state. Because the mean 
flow is baroclinically unstable and the mountain is a continous source of 
disturbances, normal-mode growth is anticipated. Therefore, a strictly 
time-independent mountain mean state is not achieveable. However, a 
"locally" steady state can be derived by limiting the analysis to that 
part of the solution in which the normal-mode amplitudes are negligibly 
small. 

In this study, the locally steady, mountain mean state is derived 
from a model run in which the Control mean flow, Y, is integrated in the 


presence of the mountain. The mountain in this case is raised from zero 
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to its full height in 12 hours. The solution is then time-averaged over 
12 hours (using fields from t = 18 to t = 30h) to filter out inertia- 
gravity waves, and this time-averaged solution is defined as the mountain 
mean state, YM. 

Test integration of Ym revealed the presence of a weak transient 
solution and some amplification of the forced modes at t = 24 h. These 
effects are present in the Interaction run which includes fm as part of 
the initial conditions for a model integration. They are included in the 
Non-interaction run by integrating Ym separately in the presence of the 
mountain and using this time-dependent Yy(t) in (6.4). Formally, then, 


the Non-interaction solution is written 


(626) 


3: Description of Experiments 


The numerical experiments conducted in this portion of the study 
are listed in Table 6.1. In each, an Interaction run and a Non-inter- 
action run are generated. Two basic states are considered: one in which 
the jetstream halfwidth, y (in (4.1) - (4.5)), is 169 latitude and one in 
which y = 80. Vertical cross-sections of the mean wind for these two 
basic states are shown in Figure 4.1. As indicated earlier, the 89 jet 
more closely resembles the upper level flow observed during cyclogenesis 
(see, for example, Palmen and Newton (1969)). However, the 89 jet also 
meets the necessary condition for barotropic instability in regions 
immediately north and south of the jet. Therefore, a 169 jet which does 


not meet this condition is included to isolate the baroclinic effects. 
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Table 6.1. List of Experiments 


Control A: No Topography 

Control B: No Topography 
I 470x240x1.5 km Mountain 
li 470x240x3.0 km Mountain 


l1 470x240x1.5 km Mountain 
IV 470x240x3.0 km Mountain 
Va 250x240x3.0 km Mountain 
Vb 570x240x3.0 km Mountain 





In these experiments, the initial amplitude of the disturbance, 
the mountain height, and the mountain width are varied. The initial dis- 
turbance amplitude ranges from very weak to approximately finite-ampli- 

tude, and is prescribed by setting tp, the time at which the perturbation 
fields are taken from the Control fun, to the appropriate value. Results 
of Trevisan (1976) suggest that lee cyclogenesis is triggered only when 
the disturbance has sufficient amplitude before encountering the 
mountain. If this is so, the results of these experiments should identi- 
fy a minimum amplitude. 

Mountain height is varied between 1.5 km and 3.0 km. The anal- 
ysis of Chung et al. (1976) indicates that mountain height plays a role 
in the generation of lee cyclones. In addition, the results of Williams 
et_al. (1981) suggest that the fundamental character of the mean flow 
distortion differs with mountain height - with air parcels tending to 
flow over smaller mountains and around larger mountains. The 1.5 km 


mountain range is representative of smaller mountains such as the 
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Appalachians, whereas the 3.0 km mountain range roughly resembles the 
Rockies. Because it represents a larger obstacle, the higher mountain is 
expected to have a greater impact on the development of the wave. 

The sensitivity of disturbance growth to mountain length is in- 
vestigated by using three mountain lengths: 470, 250, and 570. The 470 
mountain length was used in the majority of experiments and should be 


assumed unless otherwise specified. 


B- RESULTS 
1. Control Runs 

Control runs without topography were generated to: 1) measure 
growth in the absence of topography (defined here as “normal" baroclinic 
growth); 2) identify the portion of the solution dominated by linear 
dynamics; and 3) create the Y (t) fields required for initializing the 
Interaction run and defining the perturbation in the Non-interaction run. 
A control run for each basic state (y = 16% and y = 80) was required. 

Results of Control Run A (160 jet) are shown in Figure 6.1. An 
initially weak disturbance grows into a closed circulation by approxi- 
mately t = 60 h (not shown) and by t = 96 h, a mature cyclone is observ- 
ed. A log-linear plot of disturbance amplitude at wave number 8 versus 
time is given in Figure 6.2. The exponential growth predicted by linear 
theory is observed from approximately t = 30 to t = 75 h; the doubling 
time is approximately 21 hours, and is in agreement with the rough 
estimate of 24 hours for synoptic-scale baroclinic waves given by 
Haltiner and Williams (1980). Growth is not observed at other wave 
numbers until after t = 60 h, which indicates that linear processes 


dominate through this period. 
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Figure 6.la. Sea-level pressure contours for control run A at 0 h. 
Contour interval is 2.5 mb. 





Figure 6.1b. Sea-level pressure contours for control run A at 24 h. 
Contour interval is 2.5 mb. 


85 





Figure 6.1c. Sea-level pressure contours for control run A at 48 h. 
Contour interval is 2.5 mb. 





Figure 6.1d. Sea-level pressure contours for control run A at 72 h 
Contour interval is 2.5 mb. 
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Figure 6.le. Sea-level pressure contours for control run A at 96 h. 
Contour interval is 2.5 mb. 
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Figure 6.2. Amplitude of wave number 8 of v at o = 0.9167 versus time at 
450N for control run A. 
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Results of Control Run B (the 80 jet) are similar. rot in 
this case is slightly slower with a disturbance doubling time of 27 h, as 
is indicated in Figure 6.3. Even though both jets have the same maximum 
speed, the broader 160 jet leads to a larger growth rate. This is appar- 
ently because the 160 jet has more available potential energy (because of 
the larger temperature variance). Also, the stronger barotropic effects 
in the 89 jet may reduce the growth rate as has been shown by Grotjahn 


(1979). 


Time (h) 


Figure 6.3. Amplitude of wave number 8 of v at o = 0.9167 versus time at 
45%N for control run B. 
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2. Experiment I: 160 Jet, 470x240x1.5 km Mountain 


In this experiment, the effect of the 1.5 km mountain on distur- 
bance growth in the presence of the 160 jet is considered. Four cases 
are examined in which the initial wave amplitude is successively in- 
creased. In Case 1, the initial perturbation fields are taken from the 
Control at t = 24h. In Cases 2, 3, and 4, the initial fields are taken 
from 36, 48, and 60 h, respectively. In each of these cases, the distur- 
bance is initially positioned upstream of the mountain. 

In Figure 6.4, sea-level pressure contours are shown at 12-hourly 
intervals for Case 1. The disturbance crosses the mountain, and a closed 
circulation, indicative of cyclogenesis, develops on the lee side. This 
growth exceeds that observed in the corresponding control run, where a 
closed circulation is not observed until t = 60h (at t= tg + 36 h, 
relative to the time units of Case 1). As noted previously, the increas- 
ed growth is at least partially due to the superposition of the wave on 
the forced lee-side trough. To account for this afect, the non-inter- 
action run is introduced in Figure 6.4d. The closed circulation observed 
at t = tọ t 24 h is also observed here. The similarity of the two solu- 
tions suggests that the effect of the mountain on the stability of the 
mean flow is not critical to the occurrence of lee cyclogenesis in this 


Case. 


89 


Figure 6.4a. 


Figure 6.4b. 





Sea-level pressure contours for the Interaction and Non- 
interaction solutions of Experiment I, Case 1. Interaction 
solutions are given at to. Contour interval is 2/5 0mbe 
Bold contours at center of figure represent mountain. 





Sea-level pressure contours for the Interaction and Non- 
interaction solutions of Experiment I, Case 1. Interaction 
solutions are given at totl2 h. Contour interval is 

2.5 mb. Bold contours at center of figure represent 
mountain. 
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Figure 6.4c. 


Figure 6.4d. 





Sea-level pressure contours for the Interaction and Non- 
interaction solutions of Experiment I, Case 1. Interaction 
solutions are given at tot24 h. Contour interval is 

2.5 mb. Bold contours at center of figure represent 
mountain. 





Sea-level pressure contours for the Interaction and Non- 
interaction solutions of Experiment I, Case 1. Non-Inter- 
action solution is given at to+24 h. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 
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The solutions for Dase 4 are shown at tọt24 h in Figures 6.5a 
(Interaction run) and 6.5b (Non-interaction run). The Interaction 
cyclone appears more ¡intense than the Non-interaction cyclone. The 
central, closed pressure contour of the Interaction cyclone is 5 mb lower 
than in the corresponding Non-interaction cyclone. In addition, the 
pressure gradient east of the cyclone center appears stronger in the 
Interaction run. However, as in Case 1, a lee cyclone is observed in 
both runs. Similar results were observed in Cases 2 and 3 (not shown). 
Although some destabilization of the mean flow may be indicated by the 
more intense Interaction solutions, it does not appear to be critical to 
the formation of the surface cyclone. In addition, in contrast to the 
results of Trevisan (1976), the amplitude of the initial disturbance does 
not appear to be a factor in the occurrence of lee cyclogenesis. 

To isolate the effect of the mountain on disturbance growth, the 
mountain mean state, um(t), is subtracted from each of these runs to 
allow a direct comparison of yilt) and Yy(t). Recall that the latter is 
by definition simply the perturbation field from the control cm ye(t). 
A harmonic analysis (Figure 6.6) of the v field at g = 0.9167 
(p = 930 mb) along 45%N was completed for each run. The curves represent 


amplitude, A, which is computed according to 
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maigure 6.5a. 


Figure 6.5b. 





Sea-level pressure contours at t,+24 h for Interaction 
solutions of Experiment I, Case 4. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 





Sea-level pressure contours at t¿+24 h for Non-interaction 
solutions of Experiment I, Case 4. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 


93 


where a, is the amplitude of the nth harmonic. In the first three cases, 
the Interaction and Non-interaction curves are very close, which indi- 
cates that the mountain has little effect on the growth rate of the 
disturbance. In Case 4, the growth is similar through ty + 15 h, after 


which the Interaction disturbance appears to grow more slowly. 





6 12 18 


Time after to (h) 


Figure 6.6. Amplitude of wave number 8 of v at g = 0.9167 at 450N for 
Experiment I, Cases 1, 2, 3, and 4. Interaction (dashed) 
and Non-interaction (solid) solutions are shown. 


Figure 6.7, the minimum value of pg, the deviation from mean- 
state surface pressure, along 45%N is plotted versus time for each of the 
runs. As the wave moves up the mountain from t = 0 to 12 h, the Inter- 
action disturbance grows more slowly than the Non-interaction distur- 


bance. In fact, the disturbance even weakens between t - tp + 6 and 
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Figure 6.7. Minimum pg (where ps = Ps-Ps) at 450N versus time 
for Experiment I, Cases 1, 2, 3, and 4. Interaction 
(dashed) and Non-interaction (solid) solutions are shown. 


to + 12 h. As the disturbance moves down the lee side of the mountain, 
more rapid growth is observed in the Interaction run. This latter result 
Suggests that enhanced instability may exist on the lee side of the 
mountain. However, the net effect on the disturbance is smal] because of 
the period of decreased growth (or weakening) which occurs when the 
disturbance is on the windward slope. The minimum value of pg along 45°N 
for each case is listed in Table 6.2. The Interaction values are lower 


than the Non-interaction values. However, the maximum difference between 


o> 


Table 6.2. Experiment I. Cases 1-4. 
Disturbance Characteristics 


Entries in the column labeled "min pc" represent the lowest values of the 
deviation from mean-state surface pressure along 45°N at t = tg + 24h. 


Entries in the phase-speed column represent the average phase speed of 
wave number 8 of v at o = 0.9167 along 45°N over the 24-hour period of 
the model integration. 


Vm represents the orographically-forced, meridional wind component in the 
Hoskins transformation estimated from the numerical results at 
ae 


Phase differences, expressed in units of length, of wave number 8 of p< 
at 459 are given in the two columns at the right. In the subcolumn 
labeled "SG", entries represent the phase shift predicted by the semi- 
geostrophic solution (2.46) using entries from the Vm column. In the 
sub-column labeled "EXP", entries represent the differences in east-west 
position between the Interaction and Non-interaction solutions. 


Phase 
Differences 
(km) 





* 1 - Interaction; N = Non-interaction 


96 


the Interaction and Non-interaction values is less than 4 mb and does not 
appear to be critical to the formation of a closed surface pressure 
contour. 

Comparison of the phases of the two sets of solutions reveals 
that the Non-interaction cyclones lag the corresponding Interaction 
cyclones. This effect is seen more clearly by comparing the phases of 
the two solutions. In Figure 6.8, an acceleration of the Interaction 
disturbance is apparent even at t = tp + 6h, and by t = tg + 24 h, the 
phase difference between the Interaction and Non-interaction solutions 
has increased significantly. A similar phase difference is observed at 
t = ty + 24 h in Cases 2, 3, and 4 as shown in Figure 6.9. 

This result is predicted by the semi-geostrophic solution con- 
sidered in Chapter Il. It was shown there that the disturbance is 
accelerated as it passes over the mountain. Here, the Non-interaction 
run represents the linear solution (in the transform space of Chapter II) 
and the Interaction run corresponds to the semi-geostrophic solution (in 
the physical space). The phase difference between the two solutions is a 
result of the acceleration of the disturbance as it passes over the moun- 
tain. A rough estimate of this difference may be obtained using the 
transformation rule (2.46). The phase differences predicted by the semi- 
geostrophic solution and those observed (between the Interaction and Non- 
interaction solutions) in the four cases of Experiment I are given in 
Table 6.2. The value of VM required to calculate the semi-geostrophic 
value is estimated by using the maximum value of v at o = 0.9167 observed 
immediately downstream of the mountain. Qualitative agreement exists 


between the observed and theoretical values. 
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Figure 6.8 


X (x 221.27 km) 





ps (where Ps = Ps-Ps) versus X for Experiment 1, Case 1 at 
450N at 6 h intervals from tọ to tọt24 H. Dashed curves re- 
present Interaction solutions; solid curves represent Non- 
interaction solutions. East-west profile of mountain is 
shown at bottom. 
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Figure 6.9. 





AZ IT) 


ps (where Ps = Ps-Ps) versus X for Experiment I, Cases 2, 3, 
and 4 at 45%N at 6 h intervals from ty to ty+24 H. Dashed 
curves represent Interaction solutions; solid curves repre- 
sent Non-interaction solutions. East-west profile of 
mountain is shown at bottom. 
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Table 6.3. Experiment II. Cases 1-4. 
Disturbance Characteristics 


Phase 

Min ps Phase E Differences 

(mb) (m s”?) VM (km) 

Case (m s=1) 

IS N I N SG Exp 
1 -5.5 -4.2 Teed 16ed -19.6 381 175 
2 -9.1 -6.8 18.6 15.4 -23.6 459 280 
3 -14.5 all. 3 18.9 15.4 -29.6 575 290 
4 -22.6 -18.7 2 1510 -33.9 659 510 


* I - Interaction; N = Non-interaction 


3. Experiment II: 160 Jet, 479x249x3.0 km Mountain 


In this experiment, the effect of mountain height on disturbance 
growth is considered. The four cases of experiment I corresponding to 
initial disturbances from 24, 36, 48, and 60h are repeated using a 
3.0 km mountain. Results of Case 1 are shown in Figure 6.10. As in 
experiment I, the disturbance moves over the mountain and cyclogenesis is 
observed as the disturbance moves away from the lee slope. The higher 
values observed in the surface pressure ridges in this and other figures 
are due to a lack of downstream variation in the baroclinic stability of 
the mean flow and the exclusion of friction from the numerical model. In 
general, these factors will result in more intense cyclones and anti- 
cyclones in the numerical simulations, but should not change the funda- 
mental conclusions based on the comparison of the Interaction and 
Non-interaction simulations. The maximum amplitude and phase speed of 
the disturbance do not appear to be affected significantly by the higher 


mountain. The minimum value of pe in Case 1 is -5.5 mb (Table 6.3) which 
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Figure 6.10a. 


Figure 6.10b. 





Sea-level pressure contours at tg+24 h for Interaction 
solutions of Experiment II, Case 1. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 





Sea-level pressure contours at t)t24 h for Non-interaction 
solutions of Experiment II, Case 1. Contour interval is 
2.9 mb. Bold contours at center of figure represent 
mountain. 
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is close to the value of -6.4 mb observed in Case l of experiment I; the 
phase speeds are within 0.3 ms-4. In addition, the difference in 
development between the Interaction cyclone (Figure 6.10a) and the Non- 
interaction cyclone (Figure 10b) in Case 1 is less than 2 mb as in exper- 
iment I. 

The magnitude of the orographically-forced meridional wind com- 
ponent, Vm, (Table 6.3) in Case 1 is more than 6 m s-l higher than the 
corresponding value of experiment I. This increase which is also ob- 
served in Cases 2 through 4 is due to the stronger deflection of the flow 
caused by a higher mountain, and is predicted by the solution (2.13). 

Another significant difference from experiment I is observed in 
the vicinity of the mountain top at t = tọ + 24h (Figure 6.10a) in the 
form of a second cyclone. This feature first appears as a weak trough on 
the upwind slope of the mountain approximately 12 hours into the inte- 
gration and it intensifies as it moves up the mountain slope in a north- 
easterly direction. This same feature is observed in the Non-interaction 
solution (Figure 6.10b), but its amplitude is much weaker. This second- 
ary cyclone which was also observed by Walker (1982) appears to be relat- 
ed to mountain height and is discussed further in experiment Vb. 

The results of Cases 2 through 4 (not shown) are similar. In 
each case, the Interaction lee cyclone appears to be more intense than 
the corresponding Non-interaction cyclone. In Case 3, where the differ- 
ence in intensity is most noticeable, the central closed isobar of the 
Interaction cyclone is 7.5 mb lower than in the corresponding Non-inter- 
action cyclone. However, when the mean state is subtracted from the 


solutions, the difference is reduced to 3.2 mb (Table 6.3). This 
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indicates that superposition is partially responsible for the observed 
difference. In the other three cases, the minimum pg for the Non-inter- 
action run is within 4 mb of the value for the corresponding Interaction 
run. These results are in agreement with those of Experiment I, and 
again indicate that superposition is the dominant effect. 

In Figure 6.11, the amplitude of v* at % = 0.9167 along 459N is 
plotted as a function of time for each of the four cases. From t = tọ to 
t + 18H, as the disturbance moves over the mountain and down the lee 
Side, the slopes of the Interaction and Non-interaction curves are very 
close, which indicates, again, that the mountain has little effect on the 
growth rate. From t = tp + 18h to tp + 24h, after the disturbance has 
begun to move away from the mountain, some increase is observed in the 
growth rates of the Interaction disturbances. However, by this time, a 
closed circulation, indicative of cyclogenesis, has already formed and 
the additional growth does not appear to be critical to the formation of 
a closed surface-pressure contour. 

The phase acceleration noted in Experiment I is also observed 
here. The semi-geostrophic and experimental phase differences are given 
in Table 6.3. In these cases, the semi-geostrophic solution over-esti- 
mates the acceleration experienced by the Interaction disturbance. This 
is because of the large magnitude of v observed in the lee of the moun- 
tain. As a result, the semi-geostrophic solution predicts a larger phase 
acceleration. The smaller phase acceleration observed in the numerical 
Simulations may result from nonlinear blocking effects caused by the 


larger mountain. 
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6 12 18 


Time after to (h) 


Figure 6.11. Amplitude of wave number 8 of v at o = 0.9167 at 450N for 
Experiment II, Cases 1, 2, 3, and 4. Interaction (dashed) 
and Non-interaction (solid) solutions are shown. 


4. Experiment III: 80 Jet, 470x240x1.5 km Mountain 

In this experiment, the effect of the horizontal structure of the 
mean flow on the interaction between the disturbance and the mountain is 
considered. The halfwidth of the jet, Y (in 4.1 and 4.2), is reduced to 
89 and the four cases of Experiments I and II are run with a 1.5 km moun- 
tain. As noted earlier, this mean wind profile meets the necessary con- 
dition for barotropic instability both to the north and south of the jet 
Stream axis. This led to a smaller growth rate in the Control B run, and 


consequently, it is anticipated that the growth observed in this experi- 


ment will be less than that observed in Experiments I and II. 
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The Interaction and Non-interaction solutions for Cases 1 and 4 
at t = to + 24h are shown in Figures 6.12 and 6.13, respectively. As 
anticipated, the growth is somewhat less than was observed earlier. The 
maximum amplitudes of these disturbances (given as the minimum pé in 
Table 6.4) are from 3 to 15 mb weaker than those of Experiments I and II. 
In addition, a closed cyclone does not occur in every case. It is ob- 
served only in the cases in which the initial disturbance amplitude is 
large. This is in contrast to the results of Experiments I and Il in 
which the initial disturbance amplitude did not appear to be critical to 
the formation of a closed surface isobar. It is also in better agreement 
with the result of Trevisan (1976). 

As in the previous two experiments, the final amplitudes of the 
Interaction and Non-interaction disturbances are within a few millibars 
of each other in every case. However, here, the Non-interaction distur- 
bances are slightly more intense. In fact, closed circulations result in 
all but the first Non-interaction run and only in the Case 4 Interaction 
run. The amplitude of v along 450N at o = 0.9167 is plotted versul time 
for each of the four cases in Figure 6.14. During the first 15 hours, 
the growth rates of the Interaction and Non-interaction disturbances are 
in rough agreement. After t = tp + 15 h, growth in the Interaction runs 
occurs at a slower rate. 

An interesting contrast to earlier results is in the location at 
which a closed circulation is first observed. In Experiments I and II, 
the closed low was not observed until the disturbance reached the bottom 
of the mountain. Here, in Cases 2 through 4, the closed circulation is 


first observed along the lee slope close to the mountain top. 
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Figure 6.12a. 


Figure 6.12b. 





Sea-level pressure contours at t¿+24 h for Interaction 
solutions of Experiment III, Case 1. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 





Sea-level pressure contours at t¿+24 h for Non-interaction 
solutions of Experiment III, Case 1. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 
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Figure 6.13a. Sea-level pressure contours at t¿+24 h for Interation 
solutions of Experiment III, Case 4. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 
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Figure 6.13b. Sea-level pressure contours at t¿+24 h for Non-interation 
solutions of Experiment III, Case 4. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 
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Table 6.4. Experiment III. Cases 1-4. 


* I - Interaction; N 


Disturbance Characteristics 


o ese 
(m s”?) 


(m s-1) 


= Non-interaction 
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Phase 
Differences 
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Figure 6.14. Amplitude of wave number 8 of v at g = 0.9167 at 450N for 
Experiment III, Cases 1, 2, 3, and 4. Interaction (dashed) 
and Non-interaction (solid) solutions are shown. 
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Results of the phase calculations are shown in Table 6.4. As can 
be seen, the presence of the mountain causes a 20-25% increase in phase 
speed and in the Interaction runs (over the Non-interaction runs). 
However, the average values observed here are 3-5 m s-l less than were 
observed in the previous experiments. The phase lag observed here agrees 
well with the value predicted by the semi-geostrophic solution. 

5. Experiment IV: 80 Jet, 470x240x3.0 km Mountain 

In this experiment, the mountain height is raised to 3.0 km to 
examine the combined effect of the narrower jet and the higher mountain 
on disturbance growth. The first three cases of Experiment III are re- 
run for a 3.0 km mountain. In Case 4 (not shown), a large inertial 
oscillation was observed in the v fields of the Interaction run and this 
prevented a meaningful comparison with the corresponding Non-interaction 
run. 

Based on the similarity of the results of Experiments I and II, 
it is anticipated that the 3.0 km mountain will have essentially the same 
net effect on disturbance growth as the 1.5 km mountain. That is the 
final amplitudes of the Interaction disturbances should be slightly 
weaker than those of the Non-interaction disturbance and considerably 
weaker than those of Experiment II. Examination of the minimum values of 
pg for each of the three cases in Table 6.5 reveals that both of these 
expectations hold true. The Non-interaction disturbances have slightly 
lower pg than the Interaction disturbances and both sets of values are 
approximately the same as those for the corresponding runs of Experi- 
ment III (Table 6.4). In addition, the final amplitudes are from 2 to 


20 mb weaker than those of Experiment II (Table 6.3). 
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Table 6.5. Experiment IV. Cases 1-4. 


Disturbance Characteristics 


Phase 
Differences 
(km) 





x TI - Interaction; N = Non-interaction 


The Interaction and Non-interaction solutions for Cases 1 and 4 
at t = ty + 24 h are shown in Figures 6.15 and 6.16, respectively. As in 
Experiment III, a closed low is observed only in the Non-interaction runs 
of Cases 2 and 3, not in the corresponding Interaction runs. Also, as in 
Experiment III, the cyclogenesis occurs over the lee slope. 

The secondary cyclone noted in Experiment II is also observed 


here over the mountain top, although only in the Interaction runs. 
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Figure 6.154. 


Figure 6.15b. 





Sea-level pressure contours at t )+24 h for Interaction 
solutions of Experiment IV, Case 1. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. i 





Sea-level pressure contours at ty+24 h for Non-interaction 
solutions of Experiment IV, Case 1. Contour interval is 
2.9 mb. Bold contours at center of figure represent 
mountain. 
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Figure 6.16a. 


Figure 6.16b. 
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Sea-level pressure contours at t¿+24 h for Interaction 
solutions of Experiment IV, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 





Sea-level pressure contours at t¿+24 h for Non-interaction 
solutions of Experiment IV, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 
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The amplitude of the feature is less than that observed in Experiment II, 
as might have been anticipated based on the comparison of the primary 
disturbances. The lack of a secondary cyclone in the Non-interaction 
runs may be related to the proximity of the primary disturbance to the 
mountain top. In each of the Non-interaction solutions, the disturbance 
is located over the lee slopes at t = tọ + 24h. In the Interaction 
runs, it has moved over the plain well east of the mountain top. 

Results of the phase calculations are shown in Table 6.5. Phase 
speeds of the Interaction disturbances are 8-9 m s-! higher than those of 
the Non-interaction disturbances, and 5-6 m s”l higher than those of the 
Interaction runs of Experiment III. This latter result differs from that 
of Experiment II in which the increase in mountain height resulted in 
only a 1-2 m s-l increase in phase speed. This dramatic increase in 
phase speed is also evident in the large differences between the semi- 
geostrophic and experimental phase lags. It suggests that nonlinear 
effects cause additional acceleration over that predicted by the semi- 
geostrophic solution. 

6. Experiment V: Short Versus Long Mountain 

In this experiment, the effect of mountain length on disturbance 
growth is considered. The evolution of the disturbance from Case 3 of 
Experiments I and II is considered in the presence of shorter (259 lati- 
tude, north-south) and longer (579, north-south) mountains. Based on the 
results of Walker (1982), it is anticipated that less development will be 


observed with the shorter mountain than with the longer one. 


114 


The Interaction and Non-interaction solutions at t = tp + 24h 
for run Va, with the short mountain are given in Figure 6.17. Lee cyclo- 
genesis is not observed in either run. The 24-hour track of the Inter- 
action disturbance is shown in Figure 6.18. As can be seen, the distur- 
bance moves to the north of the mountain and takes a more northerly 
course than in the previous experiments where disturbances moved over the 
mountain. Lee cyclogenesis is not observed in either run. These results 
agree with those of Walker (1982) who found that disturbances tended to 
move around smaller scale mountains and exhibited little tendency toward 


cyclogenesis. 





Figure 6.17a. Sea-level pressure contours at t +24 h for Interaction 
solutions of Experiment Va, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 


LO 





Figure 6.1/b. Sea-level pressure contours at ty+24 h for Non-interaction 
solutions of Experiment Va, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure represent 


mountain. 





Figure 6.18. Position of surface-pressure trough at to, to+12 h, and 
to+24 h in Experiment Va. 
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The minimum values of pe for this experiment are given in Table 
6.6. As before, little difference in amplitude exists between the Inter- 
action and Non-interaction solutions. The secondary cyclone that was 
noted in Experiments II and IV is also observed in Figure 6.17. In the 
Interaction run, it appears as a weak low along the northern slope near 
the mountain top. In the Non-interaction run, it appears further south 
and west along the windward slope of the mountain. The phase acceler- 
ation experienced by the primary disturbance is considerable. The inter- 
action disturbance is approximately 564 km further east than the 
Non-interaction disturbance. Its phase speed is nearly 6 m gl greater 
than that of the Non-interaction run, and some 3 m s=1 greater than the 
corresponding values observed in Experiments II and IV. This rapid move- 
ment may be at least partially due to accelerations experienced by the 


Interaction disturbance along the northern slope of the mountain. 


Table 6.6. Experiment V. Cases 1-4. 
Disturbance Characteristics 


Phase 
Differences 
(km) 





* I - Interaction; N = Non-interaction 


ie, 


The Interaction. and Non-interaction solutions at t = tọ +-24 h 
are given in Figure 6.19. Lee cyclogenesis is observed in both cases, 
and the closed circulation is first observed while the disturbance is 
still over the lee slope of the mountain. This differs from the results 
of Experiments I and II in which a closed circulation was not observed 
until the disturbance reached the plain east of the mountain. The 
earlier appearance of a closed surface circulation in both the Inter- 
action and Non-interaction runs indicates that it results from the super- 
position of the disturbance on a more intense lee-side trough. This 
result, which is not predicted by the two-dimensional models considered 
in Chapter II, suggests that the north-south dimension of the mountain 
may play a role in the location and intensity of the cyclogenesis. 

Minimum pg values for this run are given in Table 6.6. As in all 
previous experiments, little difference exists between the Interaction 
and Non-interaction solutions indicating that the mountain has little net 
effect on the disturbance. Comparison of these values with those of 
Experiments II and Va yields a similar conclusion. The phase-speed cal- 
culations given in Table 6.6 show that the phase speed of the Interaction 
disturbance in this run is approximately the same as that in Experiment 
IT. The experimental phase acceleration is also much less than is pre- 
dicted by the semi-geostrophic solution. 

The secondary circulation noted in the previous 3.0 km mountain 
cases is observed here as well in the Interaction run where an intense 
small-scale cyclonic circulation exists over the mountain top. No such 
circulation is observed in the Non-interaction run, though a pressure 


trough is observed in approximately the same location. The position of 
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Figure 6.19a. 


Figure 6.19b. 





Sea-level pressure contours at t )+24 h for Interaction 
solutions of Experiment Vb, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 





Sea-level pressure contours at t )+24 h for Non-interaction 
solutions of Experiment Vb, Case 3. Contour interval is 
2.5 mb. Bold contours at center of figure represent 
mountain. 


Pug 


the secondary low-pressure center is shown at three-hourly intervals in 
Figure 6.20. The secondary cyclone is first observed at t = to + 15h as 
a weak low pressure system along the windward slope of the mountain. 
Over the next 12 hours, the low moves northeastward up the mountain, 
intensifying and broadening slightly in horizontal scale. The low 
reaches its maximum amplitude as it reaches the top of the mountain at 
t = tọ + 24 h. Thereafter, the low drifts southward along the mountain 
top and weakens. By t = tọ + 30 h, only 6 hours after it has reached its 
maximum amplitude, the secondary cyclone has disappeared. This feature 
and its behavior resemble the "type B" cyclone described by Schallert 
(1962) in his study of Colorado cyclones. It also resembles the forced 


Kelvin wave described by Gill (1977) and Murakami and Nakamura (1983). 





Figure 6.20. Position of secondary surface-pressure low at three-hourly 
intervals for Experiment Vb. 
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In this chapter, Merkine's (1975) suggestion that lee cyclogene- 
sis results from enhanced instability on the lee side of large mountain 
ranges was considered using a rough approximation of the Rocky Mountains 
as a model. Results of the numerical experiments indicate that distur- 
bances moving over mountains appear to develop more rapidly and to 
greater intensities than disturbances moving over flat terrain. However, 
these increases are shown to be due predominantly to the effect of the 
superposition of the baroclinic disturbance on the orographically forced 
lee-side trough. When the basic states are subtracted, little difference 
is observed between the amplitudes of disturbances developing on the lee- 
Side of the mountain and those developing over flat terrain. Slight 
increases in growth are observed in the Interaction runs. However, the 
effect of these increases is overwhelmed by the effect of the superposi- 
tion. This basic result is not significantly different for the experi- 
ments in which mountain height and length, initial disturbances ampli- 
tude, and the horizontal structure of the basic state are varied. In 
general, increases in the intensity of the lee-side cyclonic disturbances 
considered in this study are due to increases in either the initial 
amplitude of the disturance or the amplitude of the orographically forced 
lee-side trough, and not to enhance baroclinic instability. 

Two additional results are also noted. First, a phase accelera- 
tion is observed as the disturbance moves over the mountain. This 
acceleration is predicted by the semi-geostrophic, Eady model considered 
in Chapter Il. Qualitative agreement exists between the semi-geostrophic 
and experimental values for the 1.5 km mountain experiments, although the 


semi-geostrophic estimates tend to be larger than the experimental values 
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in the 3.0 km mountain experiments. This latter result is attributed to 
a nonlinear blocking effect caused by the larger mountain. 

Second, a secondary cyclonic circulation similar to the small-scale 
Colorado cyclone observed by Schallert (1962) and the Kelvin wave de- 
scribed by Gill (1977) and Murakami and Nakamura (1983) is observed in 
the 3.0 km mountain experiments. This feature e rapidly as it 
moves up the western slope of the mountain, but dissipates just as 
rapidly as it moves southward along the mountain ridge after reaching the 
top. At not time does this secondary cyclonic circulation appear to 


affect the development of the primary (or lee) cyclone. 
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VIT CONCLUSTONS 


The objective of this research is to explain the high incidence of 
cyclogenesis in the lee of the major mid-latitude mountain ranges. The 
study is restricted to long, north-south barriers similar to the Rocky 
Mountains. Three possible mechanisms were considered: 1) enhanced, 
lee-side baroclinic instability; 2) continuous-mode growth; and 
3) superposition. 

Simple analaytical models were used to provide a theoretical basis 
for the investigation. Eady's model was adapted to include the linear 
effects of the mountain through the lower boundary condition. The evolu- 
tion of a baroclinic wave crossing the mountain was simulated by super- 
posing the classical, time-dependent Eady solution upon the steady oro- 
graphically forced solution for Eady's model. The developing wave was 
partially cancelled as it moved up the mountain by the orographically 
forced high pressure ridge. As it reached the lee-side, however, the 
disturbance appeared to grow rapidly as it became superposed with the 
forced lee-side trough. This rapid development occurred in a manner not 
unlike that observed in lee cyclogenesis - without any actual interaction 
between the disturbance and the mountain. 

A semi-geostrophic model was considered to include finite-amplitude 
effects. It was shown that the model reduces to Eady's model if the 
Hoskins transformation is used. As with the linear model, the solution 
of this model indicated that the mountain affects the amplitude of the 


baroclinic wave only through the steady, forced solution; in other words, 
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there is no net orographic effect. However, a phase acceleration is pre- 
dicted as the disturbance crosses the mountain. For the mean flow con- 
sidered in this study, the phase acceleration, which is governed by the 
Hoskins transformation, resulted in an additional phase displacement of 
approximately 300-400 km over a 24-hour period. 

Farrell (1982) proposed that continuous spectrum growth could be 
important for lee cyclogenesis. He showed that rapid initial growth is 
possible with the appropriate initial conditions, but that the long term 
behavior gives dampening proportional to 1/t% (Pedlosky, 1964). In this 
thesis, the initial disturbance was chosen to cancel the forced mountain 
solution. The numerical experiments with the linearized spectral model 
showed that this initial disturbance was oscillatory in time, and there 
was no indication of any tendency for amplification. In experiments with 
the fully nonlinear, potential enstrophy conserving UCLA model, distur- 
bances forced by the mountain did not grow until some 30-40 hours after 
initial excitation, at which time discrete-mode growth began to dominate. 

Merkine (1975) proposed that lee cyclogenesis is more frequent 
because the lee sides of. major mountain ranges are regions of enhanced 
baroclinic instability. In this study, the fully nonlinear UCLA model 
was used with a vertically and horizontally sheared westerly flow across 
a long north-south barrier to determine the relative importance of the 
enhanced baroclinic instability mechanism and the superposition mechan- 
ism. Initial disturbance amplitude, mountain height and length, and the 
T structure of the mean flow were varied. In each experiment, 


three simulations were generated: a Control run (without topography), an 
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Interaction run, and a Non-interaction run. In the Interaction runs, the 
perturbation fields from the Control run were superposed on a numerically 
derived mountain mean state, and these conditions were integrated forward 
in time without any imposed dynamical constraints. In the Non- 
interaction runs, a time evolution of the same disturbance was 
constructed without model integration by superposing upon the Mountain 
mean state perturbation fields taken from the Control run at various 
Stages of development. These three runs were intercompared to estimate 
the growth due to normal baroclinic processes and that due to the 
Superposition of the disturbance on the orographically forced solution. 
In this way, any additional growth resulting from destabilization of the 
mean flow could be isolated. 

The results indicate that the distortions of the mean flow caused by 
the mountains do not result in significantly enhanced baroclinic insta- 
bility for the cases studied. Disturbances in the distorted mean flow 
developed at approximately the same rate as those in a parallel mean 
flow. In addition, the amplitude of the disturbance as it begins to move 
away from the mountain was approximately the same in both the Interaction 
and Non-interaction runs. This basic result is not significantly 
affected by varying the initial disturbance amplitude, the mountain 
height or length, and the horizontal structure of the mean flow. Conse- 
quently, superposition dominates over enhanced baroclinic instability in 
this study. 

Other features of the numerical solutions are also of interest. The 
mountain-induced divergent flow caused an acceleration of the disturbance 


as it moved over the mountain. This effect was predicted by the solution 
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of the semi-geostrophic, Eady model and was observed in the Interaction 
simulations of Chapter VI. The accelerations caused by the lower moun- 
tains agreed well with those predicted analytically; however, those 
caused by the higher mountain were overestimated by the semi-geostrophic 
model. This latter result is an indication that nonlinear effects are 
aa and may be a consequence of a tendency for air to flow around 
rather than over the larger barrier. 

A secondary cyclonic circulation which is not predicted by the ana- 
lytical models was also observed in the simulations of Chapter VI. This 
feature which is similar to the forced Kelvin wave described by Gill 
(1977) and Murakami and Nakamura (1982) appeared in all simulations with 
the higher mountains and in none of those with the lower mountains. The 
secondary cyclone developed rapidly as it moved up the western slope of 
the mountain and weakened, just as rapidly, as it drifted southward after 
reaching the mountain top. At no time did this feature appear to affect 
the development of the primary cyclonic disturbance nor did it show any 
tendency to develop into a lee cyclone. However, the rapid development 
of this forced disturbance suggests that enhanced instability may exist 
over the mountain for scales shorter than the synoptic scale. Its de- 
pendence on mountain height suggests that the feature is an orographic- 
ally forced circulation which depends critically on mountain size. 

Mountain length had an effect on disturbance intensity and rate of 
growth. The longer mountain, because it represented a larger barrier to 
the westerly mean flow, resulted in a larger amplitude forced ridge- 
trough pattern. This, in turn, affected the observed amplitude and rate 


of growth of the lee-side disturbance. Cyclogenesis was not observed in 
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the experiments with the shortest mountain indicating that the mountain 
length may play a role in whether or not a closed circulation is 
observed. 

For the two mean flows considered, the horizontal structure of the 
mean flow appeared to affect the growth rate of the disturbance. The 
narrower jet resulted in a lower growth rate in the Control cases of 
Chapter VI. This effect was also observed in the experiments with moun- 
tains and appears to result from there being less available potential 
energy in the narrower jet and possibly from barotropic effects. 

This research has shown that superposition is the dominant mechanism 
of lee cyclogenesis for long mountain ranges which are perpendicular to 
the basic current. When the time mean flow is subtracted, the resulting 
disturbances show little effect as they move over the mountain barrier. 
However, when the total flow is viewed, rapid cyclogenesis occurs as the 
disturbance moves down the lee-side slope under the influence of the 
mountain-forced horizontal convergence. The superposition mechanism can 
be expected to affect cyclogenesis whenever there are synoptic-scale 
forced fields. In particular, the heating off the east coasts off mid- 
latitude continents could produce a surface low which would increase the 
frequency of cyclogenesis in those regions. 

In addition, this research has shown that the amplitude of the oro- 
graphically-forced, lee-side trough is a significant factor in the 
occurrence of apparent cyclogenesis in the lee of a long meridional 
barrier. It has further shown that the amplitude of the forced trough is 


affected by the height and length of the mountain as well as the strength 
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of the current perpendicular to the mountain. These results have appli- 
cation to operational forecasting of lee cyclogenesis because they 
suggest a more important role is played by the lee-side trough. 

Hovanec and Horn (1975) have shown that static stability is decreased 
in the lee of the Rockies just prior to cyclogenesis. Static stability 
was allowed to vary in all but the linear models of this study and large 
increases in disturbance growth were not observed. However, the effect 
of feedback through diabatic mechanisms was ignored. This effect, which 
is not thought to be significant, could be evaluated by including the 
diabatic terms in the numerical simulations. 

The procedure used in this study would be improved if true steady- 
State solutions could be found for the basic flow over the mountain 
range. Possible techniques include the method developed by Kalnay-Rivas 
(1977) and the semi-geostrophic approach used by Bannon (1984) and in 
this study. 

The research in this study should be extended to other basic flows 
and other mountain barriers. For example, Smith (1984) proposed a mech- 
anism for lee cyclogenesis in the Alps which depends on a mean flow in 
which the vertical wind shear is opposite to the surface wind component 
across the mountain range. In this case, there is a large steady state 
response and the "cyclogenesis" is not a true instability. This mechan- 


ism needs to be investigated with more realistic conditions. 
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